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ABSTRACT 


Assuming  that  a  helicopter  rotor  in  forward  flight  can  be  repre¬ 
sented  by  e  flat  circular  disc  with  an  appropriate  steady  pressure 
discontinuity  across  it,  expressions  were  developed  for  relating 
the  pressure  field  (1)  to  the  total  aerodynamic  thrust  and  the 
total  steady  pitching  and  rolling  moments  experienced  by  the  rotor 
or  (2)  to  the  time-dependent  aerodynamic  flapping  moments  at  the 
roots  of  the  blades  as  they  rotate.  The  selected  pressure  field 
(1)  satisfies  the  requirement  of  a  loss  of  lift  at  the  center  and 
periphery  of  the  disc,  (2)  accounts  for  first  and  second  azimuthal 
harmonic  variations  in  the  lift  density,  and  (3)  satisfies  Laplace’s 
equation,  which  is  the  governing  equation  for  the  linearized  steady 
flow  of  an  incompressible  fluid. 

A  digital  computing  program  was  written  for  calculating  the  velocity 
field  from  the  selected  pressure  field  and  for  calculating  the 
positions  of  streamlines  emanating  from  points  above  the  disc.  Two 
sample  cases  were  chosen,  corresponding  to  the  UH-1  rotor  at 
advance  ratios,^  ,  of  0.26  and  0.08  respectively.  The  computed 
streamlines  for  both  cases  show  the  tip-vortex  phenomenon,  indicat¬ 
ing  Prandtl-Lanchester  trailing  vortices  emanating  from  the  tip 
areas  of  the  circular  wing. 

The  geometry  of  the  vortex  wake  that  is  shed  and  trailed  from  the 
blades  was  predicted  from  the  above-mentioned  streamline  calcula¬ 
tions.  The  computed  deviations  from  a  helical  wake  were  greater 
for  the  =  0.08  case  than  for  the  /*  =0.26  case . 

In  order  to  examine  the  sensitivity  of  blade-load  calculations  to 
wake  geometry,  an  existing  computer  program  for  calculating  blade 
airloads  was  modified  to  accept  the  computed  wake  geometry,  as 
contrasted  to  the  rigid  helix  usually  assumed.  The  blade-loads 
predictions  of  this  program  using  the  computed  wake  geometry  input 
were  compared  with  the  predictions  resulting  from  a  helical  wake- 
geometry  input  and  with  experimentally  measured  airloads  for  the 
two  cases  mentioned  above.  The  comparison  indicates  that  for  the 
UH-1  two-bladed  teetering  rotor  system,  the  computed  airloads  are 
more  sensitive  to  the  wake  geometry  for  the  slower  flight  {M.  =  0.08) 
than  for  the  faster  one  0“=  °  .26).  This  weaker  dependence  at 
higher  forward  speeds  may  be  due  to  (1)  the  wake  being  blown  farther 
behind  the  rotor  and/or  (2)  the  relative  predominance  of  cyclic 
pitch  input,  tangential  velocity  variation,  and  blade  flapping  motion 
in  the  total  time-varying  environment  of  the  blades  at  a  higher 
speed . 
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The  subject  investigation  was  performed  under  the  direction  of 
Mr.  William  Nettles  of  U.  S.  Army  Aviation  Materiel  Laboratories, 
Fort  Eustis,  Virginia.  The  work  was  authorized  by  DA  Task 
1F162204A13904 . 

In  preparing  parts  of  the  report,  material  was  drawn  freely  from 
the  work  of  Mr.  R.  A.  Piziali,  listed  as  Reference  2  in  this 
report.  The  authors  are  grateful  for  the  cooperation  of  Cornell 
Aeronautical  Laboratory,  and  particularly  to  Mr.  R.  A.  Piziali 
for  making  available  the  source  decks  for  their  blade- loads 
program  BLP2. 
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CHAPTER  1 


INTRODUCTION 


1.1  A  First  View  of  the  Helicopter  Wake-  Geometry  Problem 


The  desirability  of  an  accurate  determination  of  the  wake  geometry 
for  a  helicopter  rotor  stems  primarily  from  the  necessity  of  pre¬ 
dicting  the  airloads  experienced  by  the  blades.  These  airloads,  in 
addition  to  providing  the  lift  and  propulsive  force,  are  a  major 
source  of  vibratory  excitation  and  drive  the  structurally  important 
oscillatory  blade  bending  moments.  Since  aerodynamic  forces  are 
determined  by  the  flow  field,  the  relative  flow  field  must  be  known 
if  one  is  to  satisfactorily  predict  the  airloads  on  a  prescribed 
airfoil  executing  a  prescribed  motion. 

It  has  been  demonstrated  from  fixed-wing  experience  (see,  for  ex¬ 
ample  ,  the  reduced  polar  diagrams  and  reduced  lift/ angle  of  attack 
curves  in  Reference  1),  that  for  a  finite  lifting  body,  the  velocity 
field  perturbation  resulting  from  the  distribution  of  the  "free" 
vorticity  in  the  wake  must  be  taken  into  account  in  order  to  predict 
the  aerodynamic  forces  (aerodynamic  force  is  sufficiently  definitive) 
on  the  body.  It  may  be  argued  that  this  wake- induced  velocity  is 
even  more  important  for  rotary  wings  since  (1)  rotation  causes  the 
wake  vorticity  to  remain  in  the  vicinity  of  the  (succeeding)  blades 
for  a  much  longer  time  and  (even  for  the  ’  'rhest  known  advance 
ratios)  is  not  blown  behind  the  wing  so  quickly  as  with  fixed-wing 
craft,  and  (2)  as  a  result  of  a  lower  forward  speed  in  the  case  of 
a  pure  helicopter,  the  change  in  angle  of  attack  for  a  given  down- 
wash  velocity  will  be  larger  over  most  of  the  disc  than  in  the 
fixed-wing  case. 

In  addition  to  the  determination  of  the  rotor  airloads,  the  wake- 
induced  velocity  field  is  likely  to  be  influential  in  the  determina¬ 
tion  cf  (1)  interference  effects  between  rotors,  (2)  positioning  and 
effectiveness  of  auxiliary  surfaces,  (3)  interference  with  the  fuse¬ 
lage  and  other  nonlifting  bodies,  (4)  debris  entrainment,  (5)  pro¬ 
jectile  trajectories,  etc;  hence,  the  continuing  interest  in  the 
study  of  helicopter  wake  geometry  and  in  the  resultant  wake- induced 
velocity  field. 

One  of  the  difficulties  in  an  accurate  determination  of  the  wake- 
induced  velocity  should  be  immediately  apparent:  the  interdepend¬ 
ence  of  the  wake  geometry  and  the  induced  velocity  field.  The  wake- 
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induced  velocity  field  can  be  expressed  only  as  an  integral  (over 
the  wake)  of  a  function  of  the  wake  geometry  and  the  vorticity 
strengths,  while  the  wake  geometry  is  obtained  by  time  integrations 
of  the  flow  field  along  particle  paths.  As  can  be  seen  from  subse¬ 
quent  sections,  this  one  complication  alone  is  sufficient  to  make 
the  general  wake-geometry  determination  problem  difficult.  This 
interdependence  between  wake  geometry  and  induced  velocity  field 
does  not  arise  in  the  customary  treatment  of  the  fixed-wing  problem 
because  one  assumes  there  that  the  induced  velocity  is  quite  small 
as  compared  to  the  forward  speed  and  hence  the  trailed  vorticity  is 
assumed  to  lie  on  a  "rigid,  horizontal"  surface  behind  the  wing. 

Another  complicating  factor  is  the  flexibility  of  rotor  blades. 

Most  helicopter  blades  currently  have  their  first  flapwise  natural 
frequencies  within  a  few  percent  of  rotor  speed.  Table  I  summarizes 
the  flapping  natural  frequencies  for  the  UH-1  and  H-34  rotor  systems. 
It  has  been  demonstrated  that  significant  variations  in  downwash 
exist  up  to  the  4th  azimuthal  harmonic.  Under  such  circumstances, 
the  azimuthal  variations  in  aerodynamic  loads  are  almost  certain  to 
excite  substantial  oscillatory  flapping  displacements  as  well  as 
torsional  displacements,  and  the  resulting  "plunging"  velocities 
will  affect  the  aerodynamic  loads  on  the  blades  to  the  first  order, 
while  the  torsional  displacements  (i.e.,  angle  of  attack  changes  due 
to  torsion)  will  affect  them  to  the  zeroth  order.  Thus,  the  situa¬ 
tion  is  not  one  of  a  "prescribed  blade  executing  a  prescribed  motion". 
As  a  final  note,  it  may  be  mentioned  that  the  rotor  control  settings 
and  attitude  are  themselves  affected  by  th  periodic  variations  in 
the  hub  loads . 

Figure  1  indicates  the  interrelations  of  these  various  aspects  of 
the  problem.  The  figure  can  be  looked  upon  as  an  information  flow 
diagram*,  the  rectangular  blocks  representing  various  phenomena  as 
customarily  treated  in  isolation.  Thus,  for  example,  the  airloads 
and  wake  vorticity  strengths  resulting  from  the  phenomenon  of  blade 
aerodynamics  can  be  calculated  if  the  (1)  blade  structural  displace¬ 
ments  and  motions,  (2)  control  settings,  (3)  vehicle  attitude  and 
motion,  and  (4)  wake- induced  velocity  (inflow)  are  all  known  as 
"inputs"  to  the  blade  aerodynamic  equations.  These  equations  are, 
for  unstalled  subsonic  motion,  simple  algebraic  relations.  (The 
assumption  that  makes  this  possible  is  the  customary  strip  theory 
postulation  that  each  spanwise  section  of  the  airfoil  has  a  two- 
dimensional  pattern  of  aerodynamic  behavior.)  We  are  not  so  fortu¬ 
nate  in  regard  to  some  of  the  other  phenomena;  e.g.,  the  wake- 


*By  "information",  we  mean  here  "knowledge  about  the  values  (steady 
or  otherwise)  of  specific  physical  quantities  of  interest". 
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Quantities  of  interest : 

5.  Blade  displacements 
and  motions 

6.  Vehicle  attitude  and 
motion 

7  Wake  vorticity  position 
(wake  geometry) 

Figure  1.  Schematic  Representation  of  the  Helicopter 
Rotor  Aerodynamics  Problem. 


1.  Control  settings 

2.  Airloads 

3.  Vortex  strength,  position 
due  to  unperturbed  flow 

4.  Induced  velocity  field 
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induced  velocity  relations  (i.e.,  Biot-Savart  law  integrated  over 
the  wake)  and  the  wake- transport  relations.  These  are  integral 
expressions,  as  mentioned  above. 


1.2  Some  Relevant  Physical  Phenomena 


Equations  representing  some  of  the  phenomena  shown  in  Figure  1  will 
be  stated  in  this  section.  More  specific  description  of  these  phe¬ 
nomena  as  represented  in  the  computer  program  appears  in  Chapter  3 . 


1.21  Blade  Dynamics 


The  oscillatory  response  of  a  rotor  blade  to  variations  in  airloads 
can  usually  be  satisfactorily  represented  in  terms  of  its  response 
in  each  of  a  judiciously  chosen  finite  number,  ,  of  its  normal 
modes  of  displacement  (rigid  as  well  as  elastic).  If  c\K(t)  is  the 
generalized  displacement  in  the  k*1*  normal  mode,  the  flapwise  and 
torsional  displacements  at  a  point  radius  f  can  then  be  written  as* 


N 


m 


y(r>+)  '  S 

k«1 


CD 


cikCt')  *  <9°°(r) 
K 


(2) 


where  /V)-  <9'»  is  the  mode  shape  in  the  k**  normal  mode . 


*We  are  here  assuming  the  blade  to  be  infinitely  rigid  in  the  chord- 
wise  direction.  The  computations  reported  herein  used  the  additional 
structural  simplifications  of  Reference  2,  ignoring  the  dependence 
of  torsional  displacements  on  airloads,  thus  assuming  the  torsional 
response  of  the  blades  to  be  completely  known  beforehand. 


5 


The  differential  equation  for  determining  the  generalized  dis ■ 
placements  «L(+)  is  of  the  form 


^k  ik  +  '  Ik  +  ^kuk  ‘Ik  =  Gk(t) 

•••k'ljNm 

and  ^k  are  the  generalized  mass,  structural  damping, 
natural  frequency,  and  generalized  forcing  function  respectively 
in  the  k*  mode-  The  first  three  of  these  are  dynamic  "proper¬ 
ties"  of  the  rotating  blade,  as  is  the  mode  shape,  and  can  be  pre¬ 
determined  for  each  mode.  The  forcing  function  depends  on  the 
airloads  experienced  by  the  blade: 


G^Ct)  =  ^  C  *  y<K\r^  +  ’P(r,+')*  ^ 

r, 

F  (*>*')  and  represent  the  time -varying  aerodynamic 

lift  and  pitching  moment  per  unit  span  at  radius  T  • 


1.22  Blade  Aerodynamics 


Consider  a  blade  segment  with  its  leading  edge  at  x  *  — V>  and 
trailing  edge  at  x  «  +•  b  in  a  blade-attached  coordinate 
system;  see  Figure  2(a).  Let  the  relative  airflow  have  a 
"tangential"  component  Vj  in  the  Xj  direction  and  an 
"impressed"  normal  flow  component  W(x')  in  the  z.  direction, 
w  (x)  is  the  combined  effect  of  (1)  the  component  of  V/ 
normal  to  the  blade  chord  as  a  result  of  the  geometrical  angle  of 
attack  a 4  and  from  the  local  camber  of  che  airfoil  section, 

(2)  the  component  of  the  forward  speed  normal  to  the  blade  chord 
as  a  result  of  shaft  tilt,  spanwise  blade  slope,  etc.,  (3)  in¬ 
stantaneous  flapping  velocity  of  the  blade  section,  (4)  instan¬ 
taneous  torsional  velocity  of  the  blade  section  (due  to  blade 
flexibility  as  well  as  cyclic  control  inputs),  and  (5)  induced 
velocity  due  to  the  wake  vorticity.  We  wish  to  estimate  the  aero¬ 
dynamic  forces  experienced  by  this  blade  segment. 

Since  no  through-flow  can  exist  on  the  impermeable  airfoil,  the 
presence  of  the  airfoil  must  give  rise  to  such  a  "disturbance" 
that  its  effect  will  exactly  cancel  the  "impressed"  normal 
velocity  w(Xj)  in  the  region  s  .  It  can  be  shown 

(see,  for  example,  Reference  3)  that  for  a  thin  airfoil,  this 
disturbance  can  be  represented  by  a  vortex  sheet  over  the  chord 
length,  as  shown  in  Figure  2(b).  The  distribution  of  the 
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strength 


of  this  vortex  sheet  must  be  such  that 


w  OS) 


f  rro 
1  *-Xi 


ae. 


=  o 


(5) 


(a) 


blade  element  and 
"impressed"  normal 
velocity  distribution 

w  (x^  . 


(b) 


thin  airfoil  repre¬ 
sentation  with  bound 
vorticity  distribution 


YO$. 


Figure  2  . 


Physical  and  Mathematical  Blade  Element  Characteristics. 
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In  order  to  uniquely  specify  the  circulation  around  the  airfoil 
section,  it  is  customary  to  further  require  that  the  flow  be 
smooth  over  the  sharp  trailing  edge  (Kutta  hypothesis)*.  This  can 
occur  only  if 


rot)  =  o 


(6) 


Subject  to  the  condition  (6),  equation  (5)  has  now  to  be  solved. 
The  solution,  by  Sohngen’s  formula,  is 


♦b 


(7) 


After  thus  determining  the_bound  vorticity  strength  T'Cx,')  from  the 
impressed  normal  velocity  w  (x,)  ,  it  now  remains  to  deduce  the 

pressure  difference  between  the  two  faces  of  the  airfoil.  Assuming 
that  (1)  the  disturbance  velocities  are  small  compared  to  Vj  , 

(2)  the  flow  is  irrotationa 1  outside  of  the  vortex  sheet,  and  (3) 
the  flow  is  incompressible,  the  momentum  equation  for  unsteady  flow 
can  be  reduced  to 


sr(_‘t£)+v'&(~i£)  *  ~7'  air, 


(8) 


where  p  is  the  disturbance  pressure  and  <f>  is  a  velocity 


*This  is  tantamount  to  requiring  a  finite  velocity  at  the  trailing 
edge.  It  has  been  corroborated  by  experiment,  in  that  the  pressure 
distributions  on  the  surfaces  of  the  airfoil  as  predicted  from 
such  a  choice  of  the  bound  vorticity  distribution  agree  with  experi¬ 
mental  measurements  (see,  for  example,  Reference  4). 
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potential,  such  that  the  disturbance  velocity  a  equals - 

Changing  the  order  of  differentiation  in  the  first  term  and 
integrating  from  x,«—  •»  to  a  point  on  the  upper  face  of  the 
vortex  sheet, 


-TT--*  +*•*■*« 


and  similarly  at  a  point  on  the  lower  face, 


_\j 

at  fcx, 


But  it  can  be  shown  (by  considering  an  elemental  rectangular  con¬ 
tour  around  a  point  )  that 


^  ■TxT=u“-u‘-  * 


and,  by  integrating, 


s^u-^u  =  j  y (e.,1)  it. 


Using  equations  (9)  through  (12),  one  gets  the  aerodynamic  pressure 
difference  across  the  vortex  sheet: 


PL  »  Ap  <x/>  */o[Vj  /(x,,t)+ 


which  is  the  unsteady  form  of  Zhukovskii’s  theorem. 
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The  total  circulation  T  around  the  blade  section,  the  aerodynamic 
lift  F  per  unit  span,  and  the  aerodynamic  pitching  moment  about 
mid-chord  per  unit  span  are  obtained  by  integration  over  the 

chord:  * 


V  -  \  J'CXx')  cix, 

*  t 


r 


^  Ap60dx* 


(14) 


(15) 


(16) 


The  other  "output”  from  the  blade  aerodynamics  block  in  Figure  1 
is  the  wake  vortex  strength.  The  strength  of  free  vortices 
trailed  and  shed  from  a  blade  can  be  deduced  by  logical  extrapola¬ 
tion  of  the  fixed-wing  "lifting  line”  theory.  Thus,  when  the 
instantaneous  bound  vorticity  varies  along  the  blade  span,  trailed 
vortex  strengths  are  given  by 


(17) 


where  y^Cr)  is  the  strength  of  the  trailed  vortex  sheet  and 
Tw(t)  is  the  strength  of  the  bound  vortex  at  the  spanwise  loca¬ 
tion  'C  .  Similarly,  for  a  blade  with  continuously  changing 
(temporal  change)  bound  vorticity,  the  shed  vortex  strength  is  given 
by 


(r  .O  = 


ATWi(r,1t) 

V7dt 


(18) 
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where  y8<ir,0  is  the  strength  of  the  shed  vortex  sheet  as  a 
function  of  spanwise  position  r  and  time  t  ,  and  V  is  the 
total  relative  fluid  velocity  normal  to  the  blade  leading  edge. 

When  the  strength  of  the  bound  vortex  varies  both  along  the  span 
and  with  time,  the  above  expressions  for  T't  anc^  Xs  will 
represent  the  components  of  the  strength  of  the  vortex  sheet  in 
the  chordwise  and  spanwise  directions  respectively. 


1.23  Wake -Induced  Velocity  Relations 


These  relations  are  used  for  determining  the  ’’induced”  velocity 
c^Q  at  a  field  point  Q  due  to  a  distribution  of  vorticity 
f  in  a  volume  V  •  This  is  directly  analogous  to  the  problem 
of  determining  the  magnetic  induction  field  "Sq  at  Q  due  to  a 
current  distribution  J  in  V  .  The  latter  problem  has  the 
solution  given  by  the  Biot-Savart  field  law: 


/**>  f  J«  r 
4*  Jy 


d* 


(13) 


where  V1  is  the  vector  distance  from  the  source  point  to  the  field 
point.  See,  for  example,  Reference  5. 


Analogously, 


for  our  induced  velocities, 


% 


Q 


dt 


(The  /*o  is  a  characteristic  of  the  units _used  for  the  electro¬ 
magnetic  quantities  such  that  VxB  »  AoJ  •) 


If  the  source  distribution  of  vorticity  is  over  line  segments  as 
in  our  wake  model,  our  formula  is  modified  to 


(20) 
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where  r  is  the  vector  from  cU  .  ».  to  Q  ,  and  V  is  the  scalar 
vorticity  per  unit  length  along  Si  ,  an  elemental  length  of  vortex 
filament  in  the  wake.  The  integration  is  over  L.  >  which  is  the 
aggregate  of  all  the  wake  filaments. 


1.24  Wake-Transport  Relations 


Kelvin’s  theorem  of  conservation  of  circulation  (see  Reference  6) 
states  that  for  an  idea}  fluid,  the  substantial  derivative  of  the 
circulation  rc»f  v®  around  a  closed  contour  c  is  zero. 


pi; 

Dt 


*  O 


(21) 


Now,  if  c  lies  on  the  surface  21  of  a  vortex  tube, 


r.*  £v 


At 


//(V 

s 


V 


A 

n 


as  -o 


r  22) 


using  Stoke ’s  theorem  and  the  fact  that  V  *  cj  is  at  right  angles 
to  the  unit  normal  n  on  the  surface  H  .  Therefore,  using 
Kelvin’s  theorem,  one  concludes  that  a  contour  on  a  vortex-tube 
surface  will  always  remain  on  a  vortex-tube  surface.  Using  two 
such  intersecting  surfaces,  it  is  easy  to  see  that  vortex  lines 
always  travel  with  the  fluid  (Helmholtz  theorem;  see  Reference  7) . 


We  may  consider  all  vorticity  in  the  wake  model  to  have  originated 
from  the  blades.  Thus,  if  we  can  calculate  the  positions  of  fluid 
particles  released  into  the  flo  ;  field  from  discrece  locations  on 
the  disc*,  at  all  subsequent  (for  our  purposes,  discrete)  instants, 
we  will  be  able  to  predict  the  wake  geometry  (i.e.,  time -dependent 


*The  "disc"  here  means  the  surface  described  by  the  blades  as  they 
rotate.  Due  to  the  preconing  built  into  the  blades,  the  steady 
deflection  under  load,  and  the  flapping  displacements,  this  surface 
will  not,  in  general,  be  a  plane  circular  surface. 
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positions  of  the  wake-filament  end  points  in  Figure  15(b)  ) . 

1.3  Review  of  Various  Approaches  to  Handling  Rotary-Wing 
Aerodynamics  ’  '  ’  ” 

1.31  Analytical 


Among  existing  theoretical  approaches  are  the  actuator  theories, 
the  blade  element  theories,  and  the  wake  theories. 


a.  The  early  actuator  disc  theory  for  a  propeller  assumed  a 
uniform  pressure  distribution  and  induced  velocity  over  the 
disc,  and  explained  their  interrelationship  through  a  simple 
momentum  transfer  approach,  leading  to  a  relation  like 


T  a  2irRZ/°  vi  +  vf  (23) 

where  T ,  Vj  ,  and  Vp,  are  the  total  thrust,  the  induced 
velocity  at  the  disc,  and  the  forward  velocity  respectively. 
Glauert  modified  this  assumption  of  the  uniform  pressure 
distribution  to  one  increasing  linearly  from  the  front  to  the 
rear  of  the  rotor  disc  (see  Reference  8).  A  modification  due 
to  Kinner  (Reference  9)  and  Mangier  (Reference  10)  assumes  an 
azimuthally  symmetric  lift  distribution  satisfying  the 
requirement  of  a  loss  of  lift  both  at  the  center  and  at  the 
periphery  of  the  disc. 

All  the  actuator  disc  theories  implicitly  assume  an  infinite 
number  of  blades.  The  inherently  unsteady  nature  of  the 
aerodynamics  associated  with  a  finite  number  of  blades  is, 
therefore,  ignored,  as  are  all  effects  concerned  with  the 
shape  of  airfoil  sections. 

b.  The  blade  element  theories,  on  the  other  hand,  consider 
the  forces  experienced  by  individual  blades  in  their  motion 
through  the  fluid;  they  are  thus  intimately  concerned  with 
airfoil  section  characteristics.  The  blade  is  divided  into  a 
number  of  radial  elements,  and  airfoil  properties  and  momentum 
balance  are  used  for  each  element  to  yield  the  thrust  and 
induced  velocity  (see  Reference  11).  An  interesting  combina¬ 
tion  of  momentum  and  blade-element  theories  for  forward  flight 
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has  been  suggested  by  Wood  and  Hermes  (Reference  12),  They 
postulate  that  the  instantaneous  induced  velocity  at  a  point 
on  the  disc  can  be  obtained  by  superposition  of  transient 
momentum  theory  induced  velocity  distributions  corresponding 
to  the  number  of  blade  passages  over  that  point  in  space. 

c.  In  contrast  to  the  above,  the  early  vortex  wake  theories 
represented  the  wake  by  a  series  of  circular  vorticesTT  The 
geometry,  however,  was  based  on  an  assumed  induced  velocity 
field  and  on  an  infinite  number  of  blades.  The  wake  vortex 
theories  then  attempted  to  obtain  the  induced  velocities  with 
the  aid  of  the  Biot-Savart  formula.  Recent  modifications 
of  the  vortex  wake  theories  assume  helical  vortices  trailed 
from  near  the  blade  tips  and  a  "warped -pie”  type  of  mesh- 
vortex  system  that  is  shed  from  the  blade.  These  vortex  wake 
theories  may  be  looked  upon  as  a  modification  of  Prandtl's 
theory  for  a  finite  fixed  wing. 


1 . 32  Expe  rimental 


Considerable  experimentation  has  resulted  in  a  good  deal  of  test 
data  both  for  structural  loads  and  aerodynamics  of  rotary  wings. 
Excellent  smoke-trail  photographs  for  a  model  rotor  in  a  wind 
tunnel  were  obtained  by  Tararine  (Reference  13).  Very  extensive 
blade -load  and  blade -motion  data  were  presented  by  Scheimann 
(Reference  14)  and  by  Burpo  and  others  (Reference  15).  Much  of 
the  experimentation  has  been  carried  out  primarily  to  gain  insight, 
and  rigid  model  rotors  have  frequently  been  used  (e.g.,  Tararine). 

A  recent  water-tunnel  study  of  rotary-wing  flow  visualization  was 
carried  out  by  Lehman  (Reference  16).  An  interesting  electromag¬ 
netic  analog  was  reported  by  Gray  (Reference  17). 


1.33  Computational 


There  appear  to  be  very  few  numerical  solutions  of  the  simultaneous 
equations  of  wake  motion  resulting  from  blade  loading,  vortex 
strength,  and  vortex- induced  velocity  as  described  by  the  Biot-Savart 
law.  One  such  calculation  was  accomplished  by  Crimi  (Reference  18) 
with  an  extensive  iterative  calculation  for  the  time-varying  flow 
and  the  vortex  geometry. 

Reference  2  (and  its  subsequent  extension,  Reference  19)  is  an 
example  of  the  elaborate  blade-loads  programs  that  have  been 
developed  in  recent  years.  Reference  19  numerically  simulates 
most  of  the  phenomena  shown  in  Figure  1.  (A  notable  exclusion  is 
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the  wake -geometry  computation.)  We  will  (’escribe  some  details  of 
the  Reference  2  program  in  Chapter  3  of  this  report.  Another 
example  is  the  program  reported  in  Reference  20,  which  compared 
the  blade-load  predictions  using  the  rigid  wake-induced  velocity 
with  similar  predictions  using  a  uniform  inflow  velocity  at  the 
disc  for  a  four-biaded,  fully  articulated  rotor. 


1.34  Some  Complicating  Factors*  Encountered  in  the  Approaches 
“TTitecl 


1.341  Physical  Aspects 


As  a  result  of  the  inherently  unsteady  aerodynamics  of  the  rotor 
in  forward  flight,  the  "shed”  vortices  are  likely  to  be  signifi¬ 
cant,  and  this  requires  that  considerable  detail  be  retained  in 
the  wake  model  (see  Figure  15(b)  ).  Also,  due  to  rotation  the 
wake  geometry  is  more  complex,  than  for  a  fixed  wing  with  a  similarly 
varying  bound  vorticity. 

The  flexibility  of  the  blades,  as  mentioned  earlier,  makes  unpre¬ 
scribed  blade  motions  inevitable. 

"...For  example,  a  typical  blade  of  a  flapping,  three- 
bladed,  25  foot  radius  rotor  has  a  first  bending 
frequency  of  about  2.5  times  the  rotational  speed  for 
operating  tip  speeds;  by  comparison,  a  rotating  string 
or  chain  has  a  corresponding  frequency  of  about  2.4  times 
its  rotational  speed.  Certainly,  for  airloads  with 
frequencies  of  third  harmonic  and  over,  the  blade  is 
quite  flexible  and  will  respond  elastically.  This 
interaction  between  dynamic  structural  deformations  and 
airloads  is  the  second  serious  physical  complication 
in  the  blade  airloads  problem."  (Reference  21) 

1.342  Theoretical  Difficulties 


Even  if  the  equations  were  not  coupled  as  regards  blade  motion  and 


*Material  for  this  section  has  been  drawn  from  Reference  21. 
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ae rodynamics ,  the  problem  is  difficult  to  solve  because  of  the 
interdependence  of  the  wake  geometry  and  the  induced  velocity. 

To  avoid  this  complication,  drastic  assumptions  have  been  made 
regarding  vortex  geometry.  Momentum  theories,  which  bypass  wake 
effects  completely,  have  been  subject  to  other  limitations.  The 
Kinner-Mangler  theory,  for  example,  which  is  one  of  the  more 
sophisticated  actuator  disc  theories,  ignores  the  nonuniformity  of 
lift  with  azimuthal  position. 

Another  disturbing  factor  is  the  questionable  validity  of  small 
perturbation  assumptions  in  linearizing  the  basic  flow  equations. 

For  such  assumptions  to  be  valid,  the  downward  velocity  induced 
by  the  lifting  surface,  which  from  the  momentum  point  of  view  is 
the  basic  source  of  lift,  must  be  very  small  compared  with  forward 
speed.  This  is  the  case  for  a  fixed  wing  and  makes  possible  the 
assumption  that  the  trailed  vorticity  lies  on  a  rigid  cylindrical 
surface  through  the  trailing  edge  with  generators  parallel  to  the 
direction  of  flow.  It  is  much  less  obvious  that  this  assumption 
will  give  good  results  when  using  an  "actuator  disc"  to  represent 
a  lifting  rotor  at  any  forward  speed.  Although  the  downward  flow 
component  near  the  blade  tip  may  be  small  compared  to  the  total 
forward  velocity  of  the  tip,  there  are  likely  to  be  points  on  the 
rotor  where  this  is  not  the  case.  In  any  event,  the  "actuator  disc" 
car. not  recojnize  blade  rotation,  but  only  free-stream  velocities, 
so  that  "perturbations"  may  not  be  small,  and  linearized  flow 
rquat  ions  could  lead  to  large  erx>rs. 


I  .  livpr  t imental  and  Computational  Difficulties 


An  >  •  |  ii'  i  l  men'  i 1  study  of  blade  airloads  has  to  be  very  elaborate, 

In  Mm'  (ii  It  is  difficult  to  design  a  blade  which  is  both  aero- 
'yuiMij  .illy  • « r 1. 1  dynamically  similar  to  the  full  scale  item; 
i  i •  •  •  »•  r i •  |  ••  measurements  have  to  be  made  to  define  streamlines, 

in.  I*  inn,  oscillatory  blade  pressures,  drag,  and  power; 

;  '  i in i 1 1. 1 1  ton  ul  tunnel  wall  effect  is  especially  difficult  where 

.  I,,  i  I  s' r  Ilmt  ion  over  the  rotor  disc  and  unsteady  effects 

j  :J  1 1  ,  •  I  t  ,  1 1 1 1  . 


A  wi  I  (cimrt  ly  computation  such  as  that  done  in  Reference  18, 

,i  Ml,, si  ii,  ic, i'  i  1.  if  .  Is  lengthy  even  or  a  high-speed  computer  because 
,,t  i  he  1  1  me  integral  equations  (involving  unprescribed  velocities) 


I  .  4  Mot  i  .'-it  ion  ior  the  Present  Work 


The  'moke-trail  photographs  taken  by  Tararine  (Reference  13)  for  a 
riqi  :  model  rotor  in  a  wind  tunnel  showed  significant  departures 
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from  the  traditionally  assumed  helical  tip  trail,  especially  over 
the  rear  of  the  disc.  Numerical  results  of  blade-loads  computations 
(at  the  Cornell  Aeronautical  Laboratory*)  reported  by  Piziali  in 
Reference  2  suggested  that  a  more  accurate  description  of  the  wake 
geometry  might  lead  to  a  better  correlation  between  the  computed 
and  measured  blade  loads.  A  further  modification  of  the  Reference  2 
work  was  performed  by  Chang  at  CAL  (Reference  19),  who  incorporated 
(1)  the  adjustments  in  the  pitch  control  settings  and  shaft  tilt 
for  the  given  flight  condition  (thus  "representing”  the  pilot  in 
the  steady-state  computational  algorithm)  and  (2)  torsional  degrees 
of  flexibility  of  the  blades  into  the  computer  program.  Both 
Reference  2  and  Reference  19  used  a  rigid  helical  representation 
for  the  wake,  however.  The  results  of  Reference  19 

"...are  not  significantly  closer  to  the  measured  results 
than  those  obtained  in  Reference  1.**  Therefore,  it  appr  ' 
that,  at  least  for  these  configurations,  further  improve.  ..ic 
in  the  correlations  will  probably  require  the  establishment 
of  a  better  wake  model  and  an  even  more  complete  representa¬ 
tion  of  the  blade  motions...."  (Reference  19) 

Further  evidence  of  the  possible  importance  of  wake  geometry  is 
reported  in  Reference  20.  This  comparison  of  the  results  of  two 
blade -load  computations,  one  using  a  uniform  inflow  and  the  other 
using  a  variable  inflow  resulting  from  a  CAL  type  rigid  wake,  showed 
marked  differences  in  the  two  blade-load  predictions. 

It  was  concluded  that  a  blade-load  computation  using  a  more  realis¬ 
tic  wake-geometry  input  for  calculating  induced  velocities  would  be 
worthwhile .  Since  an  accurate  wake-geometry  computation  (such  as 
Reference  18)  from  a  time-varying  velocity  field  resulting  from 
wake  interaction,  fuselage  effects,  and  blade  bound  circulation 
would  be  quite  extensive,  a  modified  actuator  theory  approach  was 
adopted.  The  approximate  method  reported  here  uses  a  Kinner-Mangler 
type  actuator  disc  with  an  improved  lift  distribution  to  get  an 
approximation  to  the  time-averaged  induced  velocity  field.  The  lift 
distribution  is  chosen  so  as  to  realistically  represent  the  radial 
variation  typical  of  helicopter  rotors  and  to  accommodate  up  to  the 


*Comell  Aeronautical  Laboratory,  Inc.,  Buffalo,  New  York, 
hereafter  abbreviated  as  CAL. 

**Meaning  Piziali* s  work,  Reference  2,  according  to  our  list  of 
References . 
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second  harmonic  components  of  the  azimuthal  variation.  The  method 
then  predicts  the  spatial  positioning  of  individual  vortices  released 
from  a  finite  number  of  blades,  without  any  a  priori  assumptions 
regarding  the  wake  geometry. 
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CHAPTER  2 


THEORETICAL  MATERIAL 


2.1  Theoretical  Background 


Kinner  (Reference  9)  noticed  that  Laplace's  equation,  which  is 
the  governing  equation  for  the  pressure  field  in  case  of  linearized, 
steady,  incompressible  fluid  flow,  is  separable  in  the  ellipsoidal 
coordinates  and  that  it  is  possible  to  choose  the  coordinate  system 
and  the  harmonics  in  such  a  manner  that  the  solution  is  discontinuous 
across  a  finite,  flat,  circular  disc.  Thus,  he  related  the  pressure 
field  to  the  total  lift  developed  by  a  circular  wing.  He  also  gave 
approximate  values  of  the  induced  downwash  on  the  disc.  His  analy¬ 
sis  did  not  allow  any  disc-incidence  angle,  however,  nor  did  he 
calculate  the  sideways  induced  velocities.  In  addition,  the 
discontinuities  Kinner  dealt  with  yield  a  variation  over  the  disc 
for  certain  combinations  of  these  solutions  that  appears  quite  similar 
to  the  lift  distribution  on  a  typical  rotor. 

Mangier  (Reference  10)  modified  Kinner 's  work  to  include  the  disc- 
incidence  angle,  but  he  restricted  his  solutions  to  azimuthally 
symmetric  lift  distributions.  For  such  a  restricted  flight  condition, 
he  obtained  numerical  values  of  the  induced  downwash  on  the  disc. 
Figures  3  and  4  show  the  distribution  of  the  induced  velocity  over 
the  disc  for  a  disc  incidence  of  0°  and  15°  respectively.  Figure  5 
shows  the  same  results  as  Figure  4,  but  they  are  displayed  through 
the  variation  of  the  induced  velocity  along  three  longitudinal  lines 
on  the  disc. 


2.2  Capsule  Statement  of  the  Work  Done 


In  order  to  evaluate  the  importance  of  the  details  of  the  wake 
geometry,  it  was  decided  to  compare  the  blade-load  predictions  from 
an  existing  blade-load  program,  using  two  different  wake-geometry 
inputs.  The  program  used  is  that*  reported  in  Reference  2.  The  two 


*Referred  to  as  BLP2  of  CAL  (Blade-Loads  Program  2  of  Cornell 
Aeronautical  Laboratory) .  Some  details  of  this  program  appear  in 
Sections  3.1  and  3.2. 
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Figure 


0.3? 


A  Total  Thrust 

Ct  =  T’fV'AR* 

w  =  induced  velocity 
normal  to  disc 


Induced  Velocity  Over  Disc  for  Azimuthally  Symnjgtric 
Lift  Distribution,  Expressed  as  Contours  for  y  - ” . 
For  Disc  Incidence  =  0°.  (From  Mangier,  M  f 

Reference  10) 
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Figure  4 . 


a  Total  Hi  rust 

Ct  =  f"ViRa 


W  =  induced  velocity 
normal  to  disc 


Induced  Velocity  Over  Disc  For  Azimuthally  Symmetric 
Lift  Distribution,  Expressed,  as  Contours  f< or  — —  „ 
For  Disc  Incidence  =  15°.  (From  Mangier,  vfCf 
Reference  10) 
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Figure  5. 


Induced  Velocity  Distribution  According  to 
Mangier  and  Squire.  Disc  Incidence  =  15°. 
(From  Reference  8) 


wake-geometry  inputs  were  (1)  a  rigid,  helical  wake  geometry  result¬ 
ing  from  the  assumption  of  a  uniform  induced  velocity.  (This  is 
the  type  of  wake  geometry  used  in  the  CAL  computations  in  Reference 
2  and  will  be  termed  the  CAL  wake  geometry.)  (2)  A  distorted  wake 
geometry  resulting  from  the  steady  velocity  field  of  an  actuator 
disc  with  a  lift  distribution  that  includes  up  to  the  second  azimuth¬ 
al  harmonic  variations.  (This  wake  geometry  will  be  termed  the  UR 
wake  geometry.)  The  steady  velocity  field  used  here  includes  induced 
velocities  in  all  directions.  Specifically,  it  does  not  ignore  the 
sideways  velocities.  This  is  considered  to  be  a  definite  advantage, 
since  even  small  sideways  distortions  of  the  wake  may  be  of  impor¬ 
tance,  especially  where  the  trailed  vortex  from  one  blade  passes 
near  the  succeeding  blade(s). 

No  changes  were  made  in  BLP2  other  than  to  modify  it  to  suit  the  IBM 
System  360/65  at  the  Computing  Center  of  the  University  of  Rochester* 
and  to  modify  it  to  interface  with  the  UR  wake  computing  program. 

The  representation  of  rotor  dynamics,  blade  aerodynamics,  and  wake- 
induced  velocity  relations  was  retained  in  its  entirety. 

Section  2.3  describes  the  theoretical  scheme  leading  to  a  determina¬ 
tion  of  the  steady  induced  velocity  field  (used  for  calculating  the 
UR  wake  geometry)  for  a  given  flight  condition. 


2.3  Theoretical  Scheme 


Equations  of  Fluid  Motion 


Consider  a  stationary,  flat,  circular  disc,  radius  R  ,  in  an 
infinite  expanse  of  an  incompressible  fluid.  Establish  a  Cartesian 
coordinate  system  xyz  such  that  the  velocity  of  the  fluid  at  an 
infinite  distance  upstream  of  the  disc  is  —  Vw  3  at  a  disc- 
incidence  angle  ef  (see  Figure  6).  The  equations  for  incompressi¬ 
ble  fluid  flow  are 


momentum : 

continuity: 


•••  *  *\  - 
-i.j  *1,2,5 


(24) 

(25) 


*Hereafter  referred  to  as  UR. 
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Figure  6.  Cartesian  Coordinate  Systems  xyz  ("Wind  System")  and 
x'yk'  ("Disc  System")  in  Dimensional  Form. 


2 


where  a.  are  the  components  of  the  fluid  velocity,  /o  is  the 
constant  density  of  the  fluid,  and  p0  is  the  pressure.  For  steady 
flow,  the  momentum  equations  become 


1j%j  ‘-jofo,  i  -i.j- 1,2,5  (26) 

Noting  that  the  components  of  ^  are  (*Vf*  +  u,v,w)  ,  where  c^(u,v,w) 
is  the  velocity  perturbation,  and  assuming  a  lightly  loaded  rotor, 
i.e.,  assuming 


(u,v,w)  «  Vps 


(27) 


allows  linearizing  equation  (26)  to 


-  -  -jr  po,»  (28) 

Differentiating  (28)  with  respect  to 

VFSCV,),X  -  is 

But  the  left-hand  side  of  this  equation  vanishes  as  a  consequence 
of  the  continuity  condition  (25).  Hence, 


i.e. , 


°>ii 


=  o 


(29) 


which  is  the  governing  equation  for  the  pressure  field  in  a  steady 
incompressible  fluid  flow.  Hence,  the  disturbance  pressure  field 
p  must  satisfy  Laplace’s  equation. 
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2.32  Coordinate  Transformations  and  Solutions  of  Laplace’s  Equation 


Since  we  are  interested  in  obtaining  a  pressure  field  that  satisfies 
the  boundary  conditions  of  discontinuity  across  the  disc,  we  attempt 
to  transform  Laplace's  equation  in  3-space  to  a  coordinate  system 
suitable  for  that  purpose. 

Define  a  Cartesian  coordinate  system  x'  y*  z*  ,  which  is  obtained 
by  rotating  the  xy'z  -system  through  the  angle  a  about  the  y1  -axis 
(see  Figure  6).  For  convenience,  let  us  nondimens ionalize  our 
Cartesian  systems  on  the  disc  radius  R  ,  i.e.,  define  new  systems 
XYZand  X'Y'Z'  >  such  that 


(X.y'.z)  =  i-  (x.y’.O 


and 


R 

1 


(30) 


(X'XZ')  = 

The  transformations  between  the  wind  system  XYZ  and  the  disc  system 

x’y'z'  are: 


X 

y 

= 

z 

x' 

y 

- 

z' 

cos  a 

o 

-Sin  at 

o 

1 

O 

sina 

o 

cos  a 

Cos  o 

o 

sin  a 

o 

1 

O 

-Sin  c t 

O 

cost* 

x’ 

Y' 

- 

Z ' 

- 

X 

Y‘ 

Z 

(31) 


(32) 


Define  a  curvilinear  coordinate  system  ,  such  that 
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(33) 


It  may  be  notec'  that  this  coordinate  system  will  cover  the 

entire  3-space  once,  if  we  restrict  v  ,  y  ,  and  ytr  (see 
Figure  6)  to  the  ranges 


-1  £  *  ^  +  1 
O  <*  J?  oO 
Q<  ytr  <x  a-rr 


(34) 


Figure  7  shows  the  coordinate  system  viewed  in  the  X  Z. 

plane  (pitch  plane  of  the  helicopter) .  The  ■  constant  surfaces 
are  hyperboloids  of  one  sheet  and  the  y  «  constant  surfaces  are 
ellipsoids,  both  families  of  surfaces  being  azimuthally  symmetric 
about  the  Z1  axis.  yr  is  the  azimuthal  angle  measured  from 
the  negative  X'  axis,  counterclockwise  looking  along  the  plus  xf 
axis.  represents  the  two  faces  of  the  disc.  The  coordinate 

x/  changes  sign  as  one  crosses  the  disc. 


The  inverse  of  the  transformation  (33)  is _ 

V’*vf  z'  /[i-<AY'%r‘)]+ypV+Y'Vz'*;>A  4z'1 

7.-Z1  <35> 

'  V 
-  tan~ 

The  matrix  of  first  partials  of  the  Cartesian  coordinates  X  jY*,  Z 
(represented  by  )  with  respect  to  the  curvilinear  coordinates 

^  ^  (represented  by  ^  )  is 


27 


Figure  7.  Ellipsoidal  Coordinate  System,  Viewed  in  the 
X*Z#  Plane. 
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The  metric 

tensor  for  the 

*7^ 

system  is 

v#+7* 

o 

a.,  a  34m 

"J  -JjT- 

1! 

'at? 

1-/** 

O 

O 

u 

Va4-»* 

1  *7* 

O 

o 

(iVxn-y*; 

(37) 

which  shows  that  the  coordinates  form  an  orthogonal  system. 
This  makes  writing  out  Laplace's  equation  especially  easy,  since 
in  orthogonal  coordinates  (see  Reference  22), 


where  •k4-/?kk  (no  summation) . 

For  our  ellipsoidal  coordinate  system,  then,  Laplace's  equation 
vV s  o  takes  the  form 


*1r 


Letting 


O  (38) 


■  i%w  (39) 

it  is  found  that  (38)  separates  into  the  three  equations 


=  O 


(40) 


(41) 
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(42) 


n(v*  +  i) 


i 


O 


where  W  and  •n  art  the  constants  of  separation, 
valued  solutions  of  (40)  are 


The  s ingle - 


C 

3  L* 


cos 


integer 


(43) 


Equations  (41)  and  (42)  are  recognized  as  forms  of  Legendre's 
associated  differential  equation  (see  Reference  23),  and  their 
solutions  are 


gcm]  (44) 

‘  Lp"°7>.Qr<v>]  C45) 


Thus,  we  find  that  the  mos*  ^eneral  single-valued  solution  of 
Laplace's  equation  in  our  e_lipsoidal  coordinate  system  is 


q:< *> 


I  C cos  rn  yr 

On 

^  J 

^  Sin  m  yr 

where  wi  is  an  integer. 


(46) 


2.33  Choice  of  the  Pressure  Field,  Based  on  Total  Thrust  and  Net 
Moments 


We  shall  deal  with  the  disturbance  pressure  in  its  nondimensional- 
ized  form  P  ,  normalized  on  twice  the  free-stream  dynamic  pressure, 
i .e . , 


P  - 
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(47) 


The  general  solution  of  Laplace's  equation  is  given  by  equation  (46). 
Appendix  I,  which  was  prepared  using  Reference  24,  shows  some  of  the 
associated  functions  of  Legendre,  p"  and  Q*£  ,  and  their  proper¬ 

ties.  Before  concluding  that  equation  (46)  is  the  general  solution 
for  the  pressure  field,  we  must  examine  it  for  any  contradictions 
with  the  physical  restrictions  on  the  disturbance  pressure  P 

Thus,  we  notice  from  Appendix  I  that  the  would  be  infinite 

at  *>*±1  ,  since  they  contain  terms  with  denominators  that 

vanish  at  u*£1  .  This  leads  to  an  infinite  disturbance  pressure 

on  the  Z*  axis.  Hence,  we  must  discard  the  type  of 

component  from  equation  (44) ,  thus  reducing  it  to 


(48) 


This  form  of  eliminates  disturbances  of  infinite  magnitude 

along  the  rotor  shaft  axis . 

Similarly,  it  is  noted  from  Appendix  I  that  P" (  as  7— *>aO  . 

Since  this  would  mean  an  infinite  disturbance  at  infinity,  this 
condition  must  also  be  ruled  out  on  physical  grounds.  This  reduces 
equation  (45)  to 


§,.<7)  »  QjOy)  <«) 

Appendix  I  shows  that  (49)  tends  to  zero  as  y  tends  to  00  . 

Equations  (39),  (43),  (48),  and  (49)  are  combined  to  obtain  the 
general  form  of  the  disturbance  pressure: 


P  =2  eo*  "• i^J  (50) 


where  C?  and  D?  are  arbitrary  constants.  From  among  the  infinity 
of  terms  in  this  general  form,  we  will  retain  only  those  terms  that 
are  necessary  to  simulate  a  reasonable  lifting  rotor  or  disc. 
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Since  over  the  disc  area  the  lift  density  should  correspond  to  the 
difference  in  the  pressures  p  just  below  (  o,  *><0)  and  above 
(7=0,  v>o)  the  disc,  we  drop  all  combinations  of  in  (50)  that 

result  in  an  even  (*«♦«).  This  is  because  P*(y)  is  an  even  func¬ 
tion  if  (w  +  n  )  is  even. 

Let  us  set  the  practical  requirements  that  (1)  an  integration  of  the 
lift  density  over  the  disc  area  should  correspond  to  the  total 
(steady)  thrust  developed  by  the  rotor,  (2)  the  lift  density  should 
vanish  at  the  center  and  at  the  periphery  of  the  disc,  (3)  the 
expression  for  the  pressure  distribution  should  be  able  to  represent 
net  pitching  and  rolling  moments  supplied  by  the  lifting  disc  to 
the  craft,  (4)  the  pressure  distribution  should  reasonably  repre¬ 
sent  the  first  harmonic  components  of  lift  density  as  experienced 
by  lightly  loaded  rotor  blades,  and  (5)  there  should  be  some  pro¬ 
vision  for  simulating  the  second  harmonic  components  of  the  blade 
bending  moments . 


Total  Thrust 

Using  equation  (50),  the  total  thrust  T  is 


T 


*[  P  JA  -  ^  PdA 


[$>>> 


■  /°V«  2  2ir  C®  Q®  C  P®  (v)  rdr 

n  0  Uff**') 

From  equations  (33),  we  find  rl . R*n-**>  for  a  point  on  the  disc. 
Hence,  while  integrating  on  the  disc, 


rdr 

(  i—t+r-rnfftr) 
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♦J 

T  -  X!Cr,Q*('‘o)^Pn  (y)vdiD 

Noting  that  v=p*(v)  and  using  the  orthogonality  relations  in 
Appendix  I,  only  the  term  will  have  a  nonvanishing  integral, 

and  we  get 


T“  vjf5-cj 

Defining  a  coefficient  of  thrust  CT 


(51) 


r  a  1  -- 

UT  *  irpR+ja* 


(52) 


we  get 


r  -  A  Cl 
4  3** 


(53) 


n  is  the  rotor  speed  in  radians  per  second,  and  is  called 
the  tip-speed  ratio 


/*  =  2Tr 


(54) 


Thus,  the  combination  (m=0,>  M  *  1 )  is  necessary  for  getting  a  steady 
thrust  from  the  pressure  field. 

o 

From  Figure  8,  it  is  seen  that  qaJ)  does  not  vanish  at  the  center 
of  the  disc,  i.e.,  at  t)  =  l  1  .  In  order  to  have  a  nonvanishing 

thrust,  and  still  to  satisfy  condition  (2)  above,  we  must  include 
the  combination  Ot  H  s  1 )  along  with  any  other  combination(s) 
having  (m**©.)  •  Their  coefficients  will  have  to  be  so  related 
that  the  composite  lift  density  vanishes  both  at  Y'sO  and  atfeR  . 
Let  us  try  a  combination  of  and  (tm«  Ojh*3)»  Since  "lAs® Oat 

the  periphery  and  \)  *  -  1  at  the  center,  we  require  that 

cj-pf°(t1) -Q°  (io)  +  cj‘  p/  (±1)«  Qs  do)  =  O 

anc/  CyP°  (o)  *  u-jOo;  +  CvPj  to)  •  Q|  (io)  s  O 


(55) 
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on  tk*  Ai*,  /Lr-fuftf 


Figure  8 .  P?(*)  Plotted  as  a  Function  of  Radius  r  on  the 
Disc . 
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Thus,  the  combinations  (*^1  -  0)  tt  —  l)  and  ( tHsO ,  Pi  —  3)  can 
together  satisfy  conditions  (1)  and  (2)  on  page  32  if 


Rolling  and  Pitching  Moments 

Using  equation  (50),  the  rolling  moment  about  the  X  axis  is 


Using  orthogonality  of  the  circular  function  over  O  -  t.  V 

MrY>Vf*£<^(fo)*C?-j'*  p’  (VH-  Sin V  d  <(> 

n  (Lower-vffer) 

while  integrating  on  the  disc,  we  have 

^  0  tD 

r*df  |a  =^>[-r3-%)  1 1 

R3vP2sgAn  (o;  )dJ 

ti 

Noting  that  p*(*>)S  -  Jv  and  using  the  orthogonality 

relations  in  Appendix  I  only  the  (n-2)  term  will  survive  inte¬ 
gration,  and  we  ge*~ 

Mr  =-§.*1- irR3-/*V,fs  'O^  £)x'a*\s  (56) 


Defining  a  coefficient  of  rolling  moment 


Cm, 


Mr 
T*  R 


(57) 
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we  get 


5  •  * 

6  1 


Cr 


(58) 


A  similar  calculation  for  the  pitching  moment  M 

Hp  «  i  *TrR3/ov/p5  •  Cg  +)V  ax/s 

z’4  _  J. ;  cMp  'Cr 

~  "  a 1 


leads  to 

(59) 

(60) 


(61) 


Thus,  the  combination  (m*ij  n-z)  will  satisfy  condition  (3)  on 
page  32. 

Figure  9  shows  typical  plots  of  the  steady  and  first  two  harmonic 
components  of  experimentally  measured  blade  loads.  Since  the  radial 
variation  of  Pj(l))  in  Figure  8  is  quite  similar  to  the  variation 
of  the  first  harmonic  components  in  Figure  9,  it  was  decided  to 
retain  the  combination  (m-  1,o»4)  in  equation  (50).  Also  retained 
was  the  combination  (m  ■  2,  *  3)- 

Table  II  shows  the  combinations  of m  and  n  to  which  the  steady 
disturbance  pressure  field  P  was  restricted. 


2*  cas  w  m  (62) 

^7tf7 

2.34  Choice  of  the  Pressure  Field,  Based  on  Blade  Flapping  Moments 

y  0  y  0  y  J  «  ^ 

In  Section  2.33  the  coefficients  Lj  ,  L- j  ,  and  Uf  were 

calculated  from  Cyt  c*r  ,  and  CffM  •  The  latter  represent  the  "total”, 
or  integrated,  effect  or  the  discos  felt  by  the  craft,  in  the  form 
of  the  thrust,  rolling  moment, and  pitching  moment  respectively. 

While  this  sort  of  a  determination  of  the  pressure  field  may  seem 
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STEADY  COMPONENT 

1  ST  COS 
<ST  SiN 

2  WD  COS 
2  NO  SIN 


Figure  9.  Typical  Steady  and  First  Two  Harmonic  Components  of 

Experimental  Blade  Loads.  (From  Scheimann,  Reference  14 
H- 34 ,  Flight  20 ,  V  =  122  Knots.) 


TABLE  II 


THE  TERMS  RETAINED  IN  THE  EXPRESSION 
(62)  FOR  THE  DISTURBANCE  PRESSURE  P 


m 


n 


Remarks 


0  1  This  is  the  only  term  that  gives  rise  to  a 

non-zero  total  lift  (but  does  not  have  zero 
lift  density  at  the  center  of  the  disc) . 


Together  with  the  above  term,  will  give  zero 
lift  density  at  the  center  of  the  disc  if 

This  is  the  only  term  that  will  give  rise  to 
net  pitching  and  rolling  moments. 


1  4  Although  the  net  moments  over  the  disc  are 

zero  for  this  term,  the  radial  distribution  of 
pressure  due  to  this  term  is  quite  similar 
to  experimentally  observed  pressure  distri¬ 
butions  . 


2  3 


L 


To  account  for  second  harmonic  variations  in 
lift  density  and  blade- flapping  moments. 


quite  obvious,  it  certainly  is  not  the  only  way  to  relate  the 
pressure  field  to  physical  quantities .  For  a  real  helicopter  with 
a  finite  number  of  blades,  barring  the  case  for  which  experimental 
measurements  of  blade  differential  pressures  are  available  (as  in 
case  of  Reference  14  and  15),  it  may  be  quite  difficult  to  accu¬ 
rately  determine  the  steady  components  of  total  thrust,  rolling,  and 
pitching  moments,  because  of  uncertainties  in  fuselage  drag,  vehicle 
trim,  etc.  It  would  probably  be  advisable,  then,  to  be  able  to 
base  the  determination  of  the  steady  pressure  field  on  some  more 
easily  (and  confidently)  estimable  quantity.  The  flapwise  moment 
on  the  blade  could  be  one  such  quantity.  Even  for  an  experimental 
case,  it  is  probably  as  easy,  if  not  easier,  to  get  measurements 
of  the  blade- flapping  moments  as  it  is  to  measure  the  differential 
pressures.  In  any  event,  since  lift  variations  higher  than  the 
first  harmonic  do  not  result  in  any  "integrated"  effect  at  all,  it 
is  clear  that  expressing  the  terms  of  the  instantaneous 

flapping  moments  on  the  blade  would  be  helpful  in  simulating  higher 
harmonic  variations  in  lift  in  any  actuator  disc  approach.  This 
section  deals  with  this  approach. 

There  are  two  possible  ways  in  which  the  coefficients  C  £  a nd  D % 
can  be  expressed  in  terms  of  the  flapping  moments experienced 
by  the  blades  at  their  roots  as  they  rotate:  r 

(1)  Assuming  that  the  value  fafa  }p)7  of  the  theoretical 

*/»>•' 

steady  lift  density  at  a  point  fr>+)  on  the  disc  is  the  time 
average  of  the  lift  supplied  to  the  rotor  when  there  is  and  is 
not  a  blade  at  that  point,  or 

(2)  Assuming  that  the  M(?)  is  simply  the  flapping  moment  on 

a  blade,  when  it  is  at  the  azimuthal  position  ip  ,  and  subjected 
to  the  steady  lift  density/*^  Wj  at 

Sorer  Upper 

It  may  be  noted  that  with  assumption  (1),  the  total  momentum  trans¬ 
ferred  from  the  fluid  to  the  disc  will  equal  the  total  momentum 
transferred  from  the  fluid  to  the  P/l  blades.  Assumption  (2)  does 
not  satisfy  this  condition.  ° 

Making  assumption  (1),  let  us  further  assume  that  the  time  history 
of  the  blade  lift  density  at  a  pointf/;  ip)oi  the  disc  is  a  series 
of  pulses  as  shown  in  Figure  10.  77V>, ^represents  the  pressure 
difference  experienced  by  each  of  the  blades  at  /*  when  at  . 
p  (rt  jp)  represents  the  assumed  steady  lift  density  at  (Ktyr)  . 


P(r,f)  -  A  (ft  i/r )  -  fi(r, if)  *  (r,\+) 

Power  '  /owe  r 
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=  number  of  blades 
A  =  rotor  speed,  radians /sec 
b  =  blade  semichord 


assumed  pulsating  lift  density  for  a 
finite  number  of  blades 

assumed  steady  lift  density  for  an 
actuator  disc 


Fi;ure  10.  Time  Histories  of  the  Lift  Density  at  a  Point 
of  the  Disc. 
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Assumption  (1)  requires  that  the  areas  under  the  two  time-history 
curves  be  equal. 


TT(r,yr)=  £Z  .rftfrp) 


tower 


and 


R 


M(r)=  J  r.  JT (r,  y^)-  2i> -dr 

O 

“  **  fa  r*J(>  (rty) dr 

/outer  ■p 


Mb  'o 

Since  on  the  disc  r*ttR9(l- we  have  r2c/rj 1S^'-'R  1/  dl) 

/ovePfoce 

3 

M( *)  =r~’  •  pVps'f  vVFtf*  P(r,+)dT/. 

V*>  r> 


(63) 

o 


(64) 


Using  expression  (62)  f or  ,  we  get 

’fy"s  E {[c"cosm^*Jl^m*J 

*>  m,n 

•el 

•  Q*(io)’fp  "(t) }T/Vl-V*dl)J 


Evaluating  the  integrals  for  the  five  combinations  of  (m,n)  shown 
in  Ta^le  II,  this  redoes  to 


M(e)  —fy4  ■f-%c2-%cj+§i(tleu++J*'>*) 

+  ^irfdccs  2  Yejj/n  P  V)J 

o_  J  * 

Or,  remembering  that  Cj  —  ~g~  ,  we  nave 


<T°  = 


'6»* 


■57r2fv~_  •  P3 


Ma 


(65) 


rs 
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2+  Mb 

* Ho 

(66) 

Si 

■Mci 

(67) 

5  i  A/±> 

-'"si 

(68) 

'6  rr-fr/s  -*3 

A/b 

(69) 

/s7r£-f>y/3-X'* 

^C2 

A/b 

IV 

(70) 

/sv*rr£-_p>'*s 

M32 

where  I^o  »  ^ci  t^isi *  ^C£ ^ts2  are  t^ie  harmonic  components  of 


M(*)  *  Mo  +  MC1  COS  *  +  Msl  s//>+  +  Mca  cos  2  \ <r+M32  sio2^TJ) 


Equations  (65)  through  (70)  relate  the  steady  pressure  field  to  the 
flapping  moments  resulting  from  the  airloads  on  the  blades,  and  were 
so  used  in  the  UR  wake  computations  reported  herein.  These  rela¬ 
tions  w^re  obtained  under  the  assumption  (1)  of  equal  average 
momentum  transfer  on  page  39  . 

If  we  choose  the  simpler  assumption  (2)  on  page  39  ,  thus  ignoring 
the  average  momentum  transfer  requirement,  we  get 

MOr)  -f  2/>(r,yr).r-2b-dr  (72) 

o  '/over 


Since  on  the  disc  r2-P2(l-l)SJ  ,  we  have  rdr/^-R^dl)] 

.  o  "j 


-I 


louter 


M(V)  *  -/  dp(r,  tydi) 


Using  (62)  for  P  ,  we  get 


~  m  m  +1 

M (ir)=  -  bbR^p t)- Qjt-o)J -Of? fiJ) di)j 

t»,n  o'* 


Evaluating  the  five  integrals, 

M (-*■).  4t>#*:pr*yfj&(£ijrc*-fi7rc*)cos  & 

+  (£  i  iro*~  J  iVD*)s,„  -p  ~/6 cj cos  2  1r-/6D*si/)  2  yrj 


2  .i 

Noting  that  both  C«  and  C4  contribute  to  the  first  harmonic  cosine 
flapping  moment,  S  need  an  additional  ’’constraint”  to  determine 
them.  Let  us  choose  the  total  pitching  moment  Mjb  from  equation  (59) 
for  the  purpose.  Similarly,  equation  (56)  will  be  used  to  determine 
JJ2  a,nd.JJ2  .  The  results  are: 

£  4 


3 _ 

**•/’*'«•** 


”o 


(73) 


r°  *  _ 9 _  M 

V  8b-fY£-*s 


(74) 


Sc 


67r-fy£s  p3 


Si 


STT-f^'y-R 


Hr 


(75) 


(76) 
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(77) 


_ u _ 

25TT-b  Rs 


MC2 


9  n*  _  6  i _ _  M 

s  2  £5ir-bf>r*s'**  SJ 


(78) 


1 _ 

6*  b-f 


Mce 


(73) 


(80) 


Equations  (73),  (74)  and  (77)  through  (80)  do  not  satisfy  the 
average  momentum  transfer  condition. 


2.35  Induced  Velocity  Field 


We  shall  deal  with  the  induced  velocity  in  its  nondimensionalized 
form'gf^j^  w)  >  normalized  on  the  free-stream  velocity  ,  i.e., 


(81) 


Using  (81)  and  (47),  equation  (23)  becomes 


L“  1,2,3 


(82) 


which  is  the  completely  dimensionless  form  of  the  linearized  momentum 
equations .  Integrating  (82)  and  noting  that  all  disturbances 
vanish  at  an  infinite  distance  upstream  of  the  disc,  we  get  the 
induced  velocity  field 
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Q.(x,y,z)  =/  Pti($v,z)di; 


1,2,3  (83) 


y-  oo 


The  subscript  l  represents  the  Cartesian  "wind  cocrdinates"X,y, 
and  Z  .  In  order  to  evaluate  (83),  then,  for  a  field  point  \_ 
(X,y,  Z  ),  we  must  first  obtain  the  partials  jjf—  ■>  >  and  i  y 

of  the  pressure  field  P  ,  and  then  carry  out  the  integration  from 
00  to  £  2=  X  .  Using  the  chain  rule  for  first  partials, 


dp  -  i£f 

dx*~  4xl  '  J?7 


i>j,k  =1,2,3 


(84) 


where  "X.  represent  the  "wind  coordinates"  (X,y,Z  ) 


I? ^represent  the  "disc  coordinates"  (x',ytz'  ) 

represent  the  "ellipsoidal  coordinates"  (  Tl,  if  ) 


The  m.  are  available  from  equation  (32),  ^The  jfc  are  obtained 
by  inverting  the  matrix  of  first  partials/ infequation  (36): 


1ZL 


1 


fa 


Sin  V((l- t)£)(t+f*j] 


m 

cos  * 


V*+f 

-A 


(85) 


Ap 


The  ygy-  are  obtained  by  using  the  expression  (62)  for  P  and  the 
formulas  for  derivatives  of  Legendre’s  associated  functions  in 
Appendix  I: 


(86) 


m,n 
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dP 


tf[„f  £01}] 


(8 


/r>,n 


jfc  *  2m'?(T»aZ,(i'- v)[cCeetm  *-  c?in  m  *7 


(83) 


f»,n 


We  now  have  all  the  necessary  formulas  for  numerically  integrating 
equations  (83)  .  Some  simplification,  however,  is  possible  in  case 
of  the  X-component  of  equation  (83): 


(89) 


So  long  as  the  path  of  integration  does  not  pass  through  points 
where  p  is  discontinuous,  evaluation  of  the  integral  on  the  right- 
hand  side  yields 


u(x,y,z)  =  p(x,y,z )  oo) 

Since  the  only  place  where  the  horizontal  path  of  integration 
encounters  a  discontinuity  in  P  is  the  disc  surface,  (90)  holds  for 
all  points  upstream  of  the  disc.  For  a  point  ( X,y,Z  )  on  a  stream¬ 
line  which  crosses  the  disc,  we  use  the  physical  condition  that  the 
fluid  velocity  at  a  point  on  the  lower  face  of  the  disc  is  the  same 
as  that  at  its  image  on  the  upper  face  (i.e.,  the  actuator  is 
permeable)*,  and  evaluate  the  integral  for  U  in  two  parts: 

ujbjber  face  x 

u(x.y.z)  --  J4§(t.y.z)de  * /-£(*. y,z)d* 

//>  Me  troke  ^ «,  /outr  -fact 

he;. -c ,  U  (*,  y,z)  =  P(x, y, z) * Pfaftnecfon point) -P(i*fcrsect/on 

vake  “//fr  '  '  (91) 

/»ce  Scce 


*or,  "no  finite  changes  in  velocity  can  occur  in  the  infinitesimal 
time  it  takes  a  particle  to  cross  the  disc." 


or,  since  for  our  disturbance  pressure  P  of  equation  (62 ),P*-P 
we  get  W' 


U  (x,y»z) 

Poke 


P(x,yJz )-  2  P  (intersection  bo/ni) 

(outer  /oce  ' 


The  term  (intersection  point)  indicates  the  point  at  which  the 
horizontal  path  of  integration  up  to  the  field  point  (Xty,Z  ) 
crosses  the  disc. 

The  induced  velocity  components  V  and  W  are  evaluated  numerically. 

The  wake  geometry  resulting  from  this  steady  induced  velocity  field 
is  also  calculated  numerically,  as  described  in  Chapter  3. 
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CHAPTER  3 


THE  COMPUTER  PROGRAM 


3.1  General  Description 


This  section  briefly  describes  the  blade-loads  program,  BLP2  of 
CAL,  as  modified  at  UR.  This  program  used  to  calculate  blade  loads 
based  on  the  two  different  wake  geometries. 

The  BLP2  is  in  two  parts.  Part  1  first  generates  the  blade  segment 
coordinates  and  the  wake  element  coordinates  for  the  given  flight 
condition.  Using  these  coordinates.  Part  1  then  calculates  (using 
the  Biot-Savart  field  law)  the  induced  velocity  coefficients 
and  Tm  appearing  in  equations  (105)  and  (106),  Kj 

**  (mm  o,  3;  K  •!,  Nra  ;  j  - 1,  Nrm  J  /  » J,  Na  )■  These  S  and 
T  coefficients  are  then  written  in  binary  form  on  TOPE  2.  Figure  11 
shows  the  logical  flow  diagram  for  Part  1  of  the  unmodified  CAL  pro¬ 
gram. 

Part  2  uses  the  TOPE  2  data  set  generated  in  Part  1  as  input.  In 
addition,  all  the  dynamic  properties  of  the  blades  and  control 
inputs  to  the  craft  are  supplied  in  Part  2.  It  then  solves  the 
equations  of  rotor  dynamics,  blade  aerodynamics,  and  wake-induced 
effects  by  iteration.  (The  essential  elements  of  the  representation 
of  these  three  phenomena  in  the  numerical  program  are  summarized 
in  Sections  3.21,  3.22,  and  3.23  respectively.)  After  achieving 
convergence,  the  blade  loads  and  other  quantities  of  interest  (as 
well  as  .heir  harmonic  analyses)  are  written  on  TAPE  3.  Figure  12 
shows  the  major  steps  in  the  iterative  procedures  for  solving  the 
equations.  Figure  13  shows  the  flow  diagram  for  Part  2  of  BLP2.  It 
will  be  noticed  from  Figure  12  or  13  that  the  Part  2  procedure  in¬ 
volves  two  iterations,  one  ’’inside"  the  other.  The  iterative  solu¬ 
tion  of  the  simultaneous  equations  (109)  in  the  bound  circulations 
takes  place  "inside"  the  iterative  solution  of  the  blade  dynami- 
cauL  equations  (96). 

Figure  14  schematically  shows  the  BLP2  as  modified  at  UR.  By 
comparison  with  Figure  1,  it  may  be  noted  that  the  control  settings, 
vehicle  attitude,  and  motion  have  ceased  to  be  unknowns.  (This 
means  that  BLP2  ignores  the  "Pilot"  and  "Vehicle  Dynamics"  blocks 
in  Figure  1)  Similarly,  the  uniform  induced  velocity  is  needed  if 
the  CAL  rigid  wake  geometry  is  used;  and  some  information  on  the 
blade  loads  is  needed  if  the  UR  wake  geometry  is  to  be  computed. 

Since  the  modification  to  the  CAL  procedure  is  essentially  the  use 
of  an  alternative  wake  geometry,  the  necessary  changes  and/or 
additions  to  the  CAL  program  occur  in  Part  1  of  BLP2,  as  can  be 


k 
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Figure  11.  Flow  diagram  for  Part  1  of  BLP2 ,  (From  Reference  2) 
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Figure  12.  Major  Steps  in  the  Iterative  Procedure  for  Solving 

the  Equations  in  Part  2  of  BLP2.  (From  Reference  2) 
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SUBROUTINE  SUB21: 


COMPUTES  CONSTANTS  AND 
VARIABLES  USED  REPEATEDLY 
IN  THE  PROGRAM 

r . . 

SUBROUTINE ' SUB *  2 : 

COMPUTE?  «■  MATRICES, 
DIAGONAL  ELEMENTS  OF  r 
" I "  ELEMENTS  OF  T  EQUA¬ 
TION,  INITIAL  TRIAL  T's  AND 
GENERALIZED  COORDINATES 


ITM  =  1 

INDEX  FOR  THE  NUMBER 
OF  ITERATIONS  PERMITTED  ON 
THE  MOTIONS. 

MAX  =  NITM  _ 


INDICATE  T’s  FAILED  TO 
CONVERGE  AFTER  PRESCRIBED 
NO.  OF  ITERATIONS 


(SUBROUTINE). 

V  GAMMA  / 

V  a)  y 


/LID  \ 
T's 

^CONVERGE/ 


(subroutine), 

l  GMFRCE  / 

\  <»  y 


(subroutine 
V  GNCRE  i 

v  o)  y 


INDICATE  T  ’s  CONVERGED 
AFTER  ITM  ITERATIONS 


/  DID  \ 
MOTIONS 
.  CONVERGE y 


INDICATE  MOTIONS 
CONVERGED  AFTER 
ITM  ITERATIONS 


COMPUTE  INDUCED 
VELOCITIES  AT  ALL 
POINTS  OF  THE  DISC 


(SUBROUTINE). 
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Figure  13.  Flow  Diagram  for  Part  2  of  BLP2.  (From  Reference  2) 
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(See  Figure  1  for  other  quantities  of  interest.) 

Figure  14.  Schematic  of  the  Modified  BLP2. 
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noticed  from  the  lower  right-hand  portion  of  Figure  14.  Appendix  II 
shows  the  FORTRAN  IV  source  listings  for  subroutine  SUB14  of  Part  3, 
as  modified  at  UR,  and  for  subroutines  RWAKE,  COEFS,  STMLN,  VLCTS, 
CHAMP,  FRWRD,  AEVNC,  and  CRSNG,  which  comprise  the  UR  written  wake- 
geometry  program. 

Essentially  no  changes  were  made  in  Part  2  of  BLP2. 


5.2  Representation  of  Various  Physical  Phenomena  in  the  CAL 
Program 


This  section  briefly  describes  the  representation  of  the  phenomena 
of  rotor  dynamics,  blade  aerodynamics,  wake-induced  effects,  and 
wake  transport  in  the  original  BLP2.  An  alternative  to  the  CAL 
wake  transport  was  added  at  UR,  while  the  other  three  representations 
were  left  unchanged. 


3.21  Rotor  Dynamics 


For  a  numerical  procedure  in  which  all  pertinent  quantities  are 
computed  at  a  finite  number  A/a  of  discrete  instants  around  the 
azimuth  for  a  finite  number  A/r  of  blade  radial  stations ,  it  is 
convenient  to  harmonically  analyze  the  generalized  forces  G, 
mentioned  in  Section  1.21:  * 


_  ** 

Gk  (t)  =A0+  cos  +  s/o  t?X2  t)  (92) 

"  />mj  ft  ft 

k~  2,  Njn 

where  is  the  number  of  harmonics  retained.  (Ihe  maximum  number 
would  be  /V^~ • )  The  components  2 ,  are  defined 

by  equations  (4)  and  (92),  except  that  hereafter,  we  shall  neglect 
the  torsional  displacements  so  that  the  aerodynamic  pitching 

moments  will  not  do  any  "work"  on  the  blade  modes,  which  now  involve 
only  flapping  displacements.  Thus, 


(93) 


S3 


f 


where  :  ^ /y  is^he  length  of  the  j  th  spanwise  segment  about  radial 
position  r/^yj  \  is  the  flapping  displacement  at  KJ  in  the  k  th 
normal  mode  ;  < and  "Prt j  and  are  the  0  th  azimuthal  cosine  and  sine 
components  of  the  lift  density  P  at  the  j  th  s panwise  station, 
such  that 


P(rj  >t)  =  POJ+£ (P„j  COSO  fit  r  Q*j s/f?  oflt) 


(94) 


The  equations  of  motion  (3)  for  the  Nm  modes  ••’an  now  be  solved 
separately  for  each  of  the  harmonic  components  of  the  steady-state 
response,  resulting  from  the  corresponding  component  of  the  forcing 
function,  yielding: 


a 


°k 


fJL 


•  •  ■  £  =1 ,  Nm 


AnJe  Ca)i  -  n’rfl+BnJt  9J, 

(di’^V**?*  *!* 

These  define  the  generalized  flapping  displacements  a* 


(95) 


Q(*)~ao 

%  A 


"A 


cos  oflt+bfi  sir?  nflt) 


1 


...  /b  «  J ,  r/n 


(96) 


The  flapping  displacement  and  corresponding  "plunging”  velocity  of 
a  blade  :adial  station  can  now  be  computed  according  to  equation  (1) : 


Atm  ^ 

Vj(t)  yj  [a0  +Y(a„  cos  nA.t+b„  si»n/lt)J 

...  j  =  1 ,  Nr 


(97) 
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V 


(98) 


The  plunging  velocity  affects  the  quasi-steady  component  of  the 
bound  circulation,  as  can  be  seen  from  Figure  14. 


3.22  Blade  Aerodynamics 


The  blade  representation  of  Reference  2  used  trigonometric  series 
expressions  for  the  distributions  a^/StJand  T'/StJ mentioned  in  Section 
1.22: 


oo 


where  6  is  the  blade-chord  based  Glauert  coordinate  (see  Figure 
2(c))  defined  by 

x,  =  -b  eos  &  (101) 

The  subscript  tC  refers  Jto  the_  particular  discrete  disc  point  to 
which  the  ?(*)>  &■($),  tr0 ,  2^7  above  belong.  1C  ranges  from 

1  to  Mra  ,  wher e%aA^>^>.  See  Figure  15(a).  (Note  that  the 
vorticity  distribution  m  equation  (99)  identically  satisfies  the 
Kutta-condition  of  equation  (6).)  In  the  numerical  computations, 
the  infinite  series  in  equations  (99)  and  (100)  were  truncated  after 
the  third  term,  i.e.,/7  ranged  from  1  to  3.  (It  turns  cut  that  the 
higher  harmonic  components  in  ?x(j9) do  not  affect  the  theoretical 
7*  ,  p  ,  and  as  evaluated  from  equations  (14),  (15)  and  (16).) 

Substituting  equations  (99)  and  (100)  into  equatior\_(  5) ,  the 
coefficients^/^  can  be  related  to  the  components^  of  the  im¬ 
pressed  normal  velocity:  £ 
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Figure  15.  X- Subscript  Notation,  Wake  Configuration  and  Wake 
Vortex  Strengths.  (From  Reference  2) 


(102) 


The  total  bound  circulation  is  now  obtained  by  evaluating  the 
integral  in  equation  (14): 


(10  3) 


where 
section' 


is  the  measured  lift  curve  slope  for  the  airfoil 


The  numerical  program  also  limits  the  total  bound  circulation 
at  disc  point  K  to  a  maximum  value  ,  whirh  is  determined  from 
the  previously  known  stall  angle  of  attack  for  the  particular 

blade  section: 


(104) 


When  the  7V  as  computed  from  equation  (103)  is  found  to  exceed 
this  >  it  is  recognized  that  the  blade  is  stalled  at  that 
position  and  the  computed  is  replaced  by  .  The  circulatory 
lift,  as  computed  using  equation  (15),  is  augmented  at  the  stalled 
segment  by  a  cross-flow  drag  force. 

The  other  "output"  from  the  blade  aerodynamics  block  in  Figure  1 
is  the  wake  vorticity  strengths.  The  wake  configuration  used  in 
the  CAL  program  is  shown  in  Figure  15(b) .  The  wake  consists  of  a 
mesh  wake  behind  each  blade  for  A/fy  azimuthal  segments,  as  well 
as  concentrated  tip  (and  root)  vortices  behind  these  meshes.  The 
actual  locations  in  space  of  the  various  points  in  the  wake  model 
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are  determined  in  the  computations  of  the  "wake  vorticity  positions" 
see  Figure  1.  The  strengths  of  the  mesh-vortex  segments  are  related 
to  the  bound  vorticities  at  the  disc  points  as  shown  in 

Figure  15(c).  The  strength  of  a  concentrated  tip  (or  root)  vortex 
segment  is  taken  to  be  simply  the  maximum  value  of  the  bound  vortic¬ 
ity  on  the  corresponding  blade  when  at  the  appropriate  azimuthal 
position. 


3.23  Wake-Induced  Effects 


For  the  purposes  of  the  numerical  program,  the  vorticity  strengths 
on  all  the  filaments  in  the  mesh  wake  can  be  expressed  in  terms  of 
the/V^y  bound  vorticity  strengths  as  shown  above.  Further,  the 
wake-induced  normal  velocities  were  needed  only  at  the  Nra  points 
on  the  disc  in  the  form  of  components  similar  to  the  V/)^  in  equa¬ 
tion  (100).  For  a  given  wake  geometry,  therefore,  it  was  possible 
to  evaluate  coefficients  such  that  the  induced  normal  velocity 

at  the  /fth  disc  point  due  td  the  mesh  wake  is  given  by 


(e)  =  z  (s„K .  +£  .  cos  r>e)ij 

Jml  J  />*!  J 


(10  5) 

*ra 


Similarly,  for  the  concentrated  trailed  vortices  part  of  the  wake, 
coefficients  7/»^.  are  computed  such  that  the  corresponding  induced 
normal  velocity 4due  to  the  trailed  wake  is  given  by 


Mo  J  _ 

^  (0)  (T°«  ■  +L  T/iK.C0S  (106) 

£mJ  1  < 

.  . .  2,  N(r& 

Tf  is  the  maximum  bound  vorticity  on  a  blade  when  at  azimuthal 
position  i  .  The  induced  velocity  coefficients and  .  are 
evaluated  in  Subroutine  SUB14  of  Part  1  of  the  CAL  program.  1 

The  "impressed"  normal  flow  component  Utl  '(e)  is  broken  up  into  two 
parts:  one,  v«(e),  resulting  from  all  ttie  contributing  quantities 
on  page  6  except  the  .wake  vorticity;  and  the  other,  (e),  resulting 
exclusively  from  the  wake  vorticity.  Thus, 

^(s)  =  (g)  *  (e)  <  lov) 


•  •  •  ^ 1  t 
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and 


(108) 


Combining  equations  (102),  (103),  (107),  and  (108),  it  is  possible 
to  write 


£  -4 


+  Z<T«jrj 

J=1  J 


\ ra 


( i:  9) 


which  is  a  set  of  equations  for  the  NrQ  unknowns  JJ.  .  Ju  Is 
the  bound  circulation  resulting  from  the  normal  f low  ,  and 

represents  the  bound  circulation  at  disc  point  K  resulting  from  J  J 
the  inducer,  velocity  due  to  the  -wake  vorticity  segments  whose 
strengths  depend  on  the  bound  circulation  at  disc  point  j  .  It 
may  be  noted  that  the  depend  on  the  blade  response,  since  the 
V*  include  the  blade-flapping  velocities,  torsional  displacements, 
and  torsional  velocities.  Equations  (109)  are  solved  by  iteration 
in  subroutine  GAMMA  of  Part  1  of  the  CAL  program. 


3.24  Wake  Transport 


Representation  of  the  wake  transport  in  the  CAL  program  consists  of 
determining  the  positions  of  the  points  defining  the  mesh  and 
trailed  vortices  in  the  wake  model  shown  in  Figure  15.  The  rigid- 
wake  calculations  reported  in  Reference  2  used  the  following  formu¬ 
las  for  the  x7,  V  and  z'  coordinates  of  the  point  :  :  t  in  the 
wake  model,  {iieB  ):  J  <X 

*ijta  *~vrs ' SL  ' O'- f+sy  cos* - rj cosfa-i+Sj* v] 


rJ S'n[(ka-L+ cos[(4:<fl+S)A v]  (2JCN 


zijka  “  (  vrs  S''1*-#)  *  7T-  -(i  -  7 "  &)  ~yQC  (o) 


.-.I  >2 


where  is  the  azimuthal  increment,  .  Mq  is  the  azimuthal 

position  of  the  blade  shedding  the  wake.^>?  is  the  semichord  at 

the  blade  segment  end  point  j  •VoeC^'j)  is  the  built-in  coning  plus 
steady  deflection  of  the  blade  Segment  end  points,  positive  upwards. 
c£  is  the  assumed  uniform  induced  velocity  (quantity  4”  in  Figure 
14)  over  the  disc.  0  is  the  fraction  of  a  time  interval  through 
which  tfhe  wake  is  advanced  with  respect  to  the  true  blade  position 
fcci  •  w  >  which  is  quasi-empirical,  was  determined  in  such  a  way 
that  for  a  two-dimensional  oscillating  airfoil,  this  discrete 
representation  of  the  shed  wake  gave  results  that  were  in  ’’reasona¬ 
ble  agreement”  in  the  ’’reduced  frequency  range  of  interest"  with 
the  classical  solution  of  the  problem.  (See  Reference  2  for  details 
of  this  determination  of  0  .)  For  /  =  1,  the  coordinates  (xjy,r')4*;£ 
are  simply  the  points  on  the  trailing  edge  of  the  blade  at  azimuthal  a 
position  Ka  . 

It  will  be  noticed  that  this  wake  transport  ignores  the  displace¬ 
ments  in  the  sideways  direction.  Also  neglected  is  the  variation, 
over  the  disc,  in  the  induced  velocity  in  the  Z'  direction. 


3.3  Algorithm  for  Computing  the  UR  Wake  Geometry 


It^can  be  seen  from  equations  (110)  that  the  point  on  the  disc* 

Rf: j.  (*', y.z'J  where  the  fluid  particle  now  at  the  point  *ijka 

J  *  in  the  wake  used  co  be  (<■-/**•<$”)  time  intervals  before 
the  blades  assumed  the  azimuthal  position  ka  ,  has  the  coordinates 


x°/yAa  ctsfa-i+S)*  vf-  f  bj  sin[(ka-i+S)fi  +] 

_  (HD 


y°ijka  =  n  sinfa-i+f)*  *J-jrbj  COs[(ka-i+£)&  Wj 


■'■‘The  "disc”  is  not  a  flat  surface,  since  blade  steady  deflection 
and  built-in  coning  are  accounted  for,  through  ■ 
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The  difference  between  the  expressions  (110)  and  (p.11)  defines 
the  movement  of  this  fluid  particle  from  (x°>  V* *  z*')  to 

,  under  the  influence  of  Vpp  and  UT ,  in  ) 

time  intervals . 

In  the  UR  wake -geometry  computations,  the  points  of  ’’injection1’ , 

I j Jca  ,  can  be  retained*  as  in  (111),  but  the  movement  under  the 
influence  of  both  and  the  induced  velocity  was  computed  differ¬ 
ently,  through  an  actual  calculation  of  the  streamlines  through  the 
injection  points.  The  steady  pressure  field  used  for  calculating 
the  streamlines  always  had  the  pressure  discontinuity  on  the  flat 
disc . 

3 . 31  Subroutine  RWAKE 


Subroutine  RWAKE,  which  is  the  entry  point  for  calculating  the  UR 
wake  geometry,  first  chooses  the  steady  pressure  field  through 
formulas  such  as  those  developed  in  Sections  2.33  and  2.34.  Figure 
16  shows  the  flow  diagram  for  subroutine  RWAKE.  Table  III  shows 
the  various  options  for  the  choice  of  the  steady  pressure  field 
based  on  various  physical  quantities. 

It  may  be  noticed  from  Appendix  I  that  the  values  of  Q \n  V*  7)  are  m 
purely  imaginary  quantities  when  /)  is  even.  The  corresponding  £ 
and  2)™  also  turn  out  to  be  purely  imaginary  quantities.  (J®'  w 
are  always  real.)  In  such  cases,  the  quantities  actually  handled 
under  QMN  (J)  wer &  i  .  To  offset  this  factor  i  ,  the  corre-  _  m 
spondingquantities  CMN(J)  and  DMN(J)  in  core  were  set  equal  to  -i  Cn 
ar\d-i  2)f”  respectively.  Thus,  all  the  computations  were  done  in 
’’reals",  and  the  products  CMN(J)*QMN(J) ,  DMN( J)*QMN( J)  always  agreed 
with  their  theoretical  counterparts  CTQ%(ty) (tV)  f°r 
J;  i.e.,  for  all  combinations  of  ft*  ana  n. 

After  choosing  the  steady  pressure  field,  RWAKE  sequentially  calcu¬ 
lates,  using  equations  (111),  the  disc-injection  points  ^ijka  for 
the  mesh  and  trailed  wake  grid  positions  (see  Figure  15(b)  );  it 
then  calls  subroutine  STMLN ,  which  calculates  the  streamline  through 
a  given  disc  point.  At  the  end  of  RWAKE,  all  the  wake  coordinates, 
disc  injection  point  coordinates,  total  normal  velocity  values  at 
disc  points,  and  induced  velocity  values  at  the  disc  points  are 
written  on  TAPE8. 


*An  option  in  the  UR  wake  computation  program  also  enables  the 
particles  to  be  injected  on  the  flat  disc,  thus  ignoring  tyoc  • 
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TABLE  III.  VARIOUS  OPTIONS  FOR  CHOOSING  THE  STEADY  PRESSURE 
FIELD  P  THROUGH  THE  COEFFICIENTS  C ?  AND  2?T 

n 

p  -  Z  Q”(i?)fcZC0S  m  ™  #] 

_ f»»n  _ _ _ 


KKSW 


Remarks 


1  Subroutine  COEFS  bypassed.  C^and  remain 

unchanged. 


2 

^  m  fS9 

andJ/,1  to  be  read  in  as 
routine  COEFS. 

input  to  sub- 

3 

C™  and  D™  calculated  from 
A  •  (See  Section  2.33.) 

CT  ■  >  CMt ,  and 

4 

C™  and  2?^  calculated  from 
and  A  .  (See  Section  2.33) 

TT,Mj>,Mr,R,  VFS 

5 

C™  and  calculated  from 

input  in  the  form 

of  harmonic  components  of  the  flapping  moment 
experienced  by  the  blade  as  it  rotates.  (See 
Section  2.34.)  Average  momentum  transfer  con¬ 
dition  satisfied. 


6  C n  andjP„  calculated  from  input  in  the  form 

of  harmonic  components  of  the  blade  flapping 
moments,  plus  *?,/>,  V/rj  >  b,  .  Average 

momentum  transfer  condition  not  satisfied. 

7  C™  ar\dJ7%*  calculated  from  v,s>  Nb. 

plus  the  harmonic  components  of  the  blade  flap¬ 
ping  moments  which  are  computed  from  the  har¬ 
monic  components  of  the  experimentally  measured 
lift  density  on  instrumented  blades.  Average 
momentum  transfer  condition  satisfied. 


C„  and^^7  calculated  from  the  same  inputs  as 
7  above ,  plus  Mr  and  W  :  Average  momentum 
transfer  condition  not  satisfied. 


8 


3.32  Subroutine  STMLN(KEND) 


Subroutine  STMLN(KEND)  follows  the  fluid  injected  from  a  given  disc 
injection  point  for  a  number  (KEND-1)  of  time  intervals,  saving  the 
coordinates  of  the  position  after  each  time  interval.  The  field  is 
assumed  steady,  and  the  particle  paths  and  streamlines  are  identical. 
A  time  interval  is  equivalent  to  the  rotation  of  a  blade  from  one 
azimuthal  position  to  the  next.  Thus, 

At  -  seconds  (112) 

To  nondimens ionalize  time,  we  must  use  a  normalization  that  is  com¬ 
patible  with  the  normalizations  (30)  and  (81)  for  distance  and  veloc¬ 
ity  respectively.  This  calls  for  normalizing  t  on  the  time  taken 
by  Vrj  to  travel  a  distance#  ;  i.e., 


At  - 


-  &  & 


(113) 


where  Alt  is  the  angle  (radians)  between  two  azimuthal  positions  and 
/t>  is,  again,  the  tip-speed  ratio.  For  accurate  streamline  calcu¬ 
lations,  the  program  has  the  provision  for  further  reducing  the  time 
increment  through  a  factor  NINC.  Thus,  each  time  increment  TINC  is 
(l/NINC)  times  the  time  interval,  and 


T/NC  = 


(114) 


Equation  (114)  represents  the  "usual"  time  increment  value.  An 
exception  is  the  very  first  interval  of  a  streamline.  (See  Figure  17) 
For  a  fractional  wake  advance^*,  the  TINC  for  the  first  NINC  incre¬ 
ments  of  each  streamline  will  be  (/-  <f)  times  that  in  equation  (114) . 
NINC  is  one  of  the  accuracy  parameters  read  in  at  the  beginning  of 
RWAKE. 


Subroutine  VLCTS  numerically  calculates  the  induced-velocity  com¬ 
ponents  V  and  W  at  a  given  point  ( X,y,Z  )  by  the  procedure  indica¬ 
ted  in  Section  2.35.  XUP  (the  location  of  "infinity"  upstream  of 
the  disc)  and  XINC  (the  value  of  the  incremental  length  used  for  the 
numerical  integration  of  (83)  )  are  the  accuracy  parameters  used  in 
VLCTS.  Both  XUP  and  XINC  are  also  read  in  at  the  beginning  of  RWAKE. 
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*  VLCTS  calculates  V  and  W  at  any 
point  by  integration  from  an 
infinitely  upstream  distance . 

**CHAMP  calculates  P ,  ,  and -ft 

at  a  given  point.  9 


Figure  17.  Flow  Diagram  for  Subroutine  STMLN(XEND) 
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The  quantities  KFWD  and  KRVRS  appearing  in  Figure  17  are ’’counters” 
for  the  number  of  flat-disc  crossings  made  by  Che  particle  being 
followed,  in  the  forward  (i.e.,  normal,  or  downward)  and  reverse 
(i.e.,  upward)  directions  respectively. 

„  dF>  dP  , 

Subroutine  CHAMP  computes  r3  ,  £y  ,  and  TF’at  a  given  point  (V,y,V) , 

using  the  formulas  from  Section  2.35.  It  may  He  noticed  that  the 
multiplier  appearing  in  equation  (87)  for  will  make  the  quantity 
real  since,  as  explained  in  Section  3.31,  either C™ ,2™  , 

but  not  both,  will  be  imaginary  for  each  value  of  ».  The  quantity 
handled  in  QMNMl(J)  will  always  be  real,  being  equal  to  if 

n  is  odd,  and  equal  to 01?)  if  n  is  even.  his  leads  to  agree¬ 
ment  between  the  products  CMN(J) *QMNM1(J) ,  DMN ( J ) *QMNM1 ( J )  and  their 
theoretical  counterparts 


CC('7h  fo'0"T 

3.33  Subroutine  FRWRD 


Subroutine  FRWRD,  which  is  called  in  subroutine  STMLN,  computes  the 
next  position  on  the  streamline  by  invoking  the  (call  ADVNC  -  call 
CHAMP)  cycle  NINC  times.  Figure  18  shows  the  flow  chart  for  sub¬ 
routine  FRWRD.  Subroutine  ADVNC  saves  the  current  values  of 
coordinates,  velocities,  pressure,  etc.,  and  calculates  the  coordi¬ 
nates  of  the  particle  after  a  time  increment  TINC,  using  the  current 
values  of  the  velocity.  If  a  flat-disc  crossing  is  recognized  on 
ADVNC-ing,  FRWRD  compensates  for  it  by  calling  CRSNG.  Subroutine 
CRSNG  assures  that  the  values  of  the  velocity  are  continuous  across 
the  flat  disc,  although  the  pressure  has  a  discontinuity  at  the 
flat-disc  surface.  Ihis  is  achieved  in  CRSNG  by  first  finding  out 
the  exact  point  of  crossing,  and  then  modifying  the  O  component  of 
the  velocity  through  an  operation  similar  to  equation  (91). 


3.4  Integration  With  the  CAL  Program 


It  is  seen  from  Figure  11  that  the  CAL  wake  coordinates  are  generated 
in  subroutine  SUB14  of  Part  1  of  the  original  BLP2.  (SUB14  calls 
another  subroutine  COORD  for  actually  computing  these  coordinates.) 
Since  our  modifications  to  the  BLP2  essentially  consist  of  generat¬ 
ing  and  using  an  alternate  wake  geometry,  the  interfacing  between 
the  UR-written  wake -computing  part  of  the  program  and  the  BLP2 
occurs  in  subroutine  SUE14 .  The  FORTRAN  IV  source  listing  for 
SUB14  appears  in  Appendix  II.  Figure  19  shows  the  flow  diagram 
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*ADVNC  saves  the  values  of  various  quantities  and  calculates 
coordinates  of  new  position  (after  a  time  increment  of  TINC)  using 
current  velocity  values. 

**CRSNG  compensates  for  crossing  of  the  disc,  so  that  velocity 
values  on  either  side  of  the  disc  are  continuous,  although  pressure 
is  not . 


Figure  18.  Flow  Diagram  for  Subroutine  FRWRD. 
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♦KWAXE  computes  the  coordinates  of  the  UR  wake 
and  writes  them  on  1APE8 


Figure  19.  Flow  Diagram  for  Subroutine  SUB14  as 
Modified  at  UR. 
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for  subroutine  SUB14  as  modified  at  UR.  The  changes  and  addi¬ 
tions  to  the  original  program  are  shown  in  dotted  lines.  The 
changes  were  made  so  as  to  make  it  possible  to  compute  and  use 
either  of  the  two  wake  geometries.  The  UR  wake  geometry  is  gener¬ 
ated  (and  written  on  TAPE8)  in  subroutine  RWAKE.  Provision  was 
also  made  to  make  it  possible  to  generate  just  the  UR  wake  and  then 
terminate  the  job.  Conversely,  a  UR  wake  generated  in  a  previous 
run  can  be  directly  read  off  TAPE 8  and  used  for  computing  the 
induced  velocity  coefficients  in  Fart  1.  All  these  options  are 
accomplished  through  an  input  integer  variable  MMSW.  (See  the 
source  listing  for  subroutine  RWAKE  in  Appendix  II  for  a  description 
of  these  options.) 


3.5  Peripheral  Program 


Two  peripheral  programs  had  to  be  written  in  addition  to  the  UR 
wake  computation  described  in  Section  3.3.  These  were  named 
RWAKEPIC  and  RSLTSPLT.  Both  make  extensive  use  of  the  plotter 
subroutines  package  a . ailable  on  the  University  of  Rochester  disc 
library. 

RWAKEPIC  reads  the  UR-l omputed  wake  geometry  from  the  TAPE8  data 
set  generated  in  subroutine  RWAKE  of  Part  1  and,  using  the  plotting 
routines,  pictorially  displays  the  wake  geometry.  Three  options 
are  available  as  to  the  type  of  picture  produced:  pictures  of 
streamlines,  instantaneous  pictures  of  the  wake  generated  by  a  fi¬ 
nite  number  of  blades,  and  progressive  distortion  of  stream- tubes 
originating  (from  points  equispaced  around  the  azimuth)  on  the  disc. 
RWAKEPIC,  which  is  a  separate  main  program,  was  found  to  be  useful 
in  interpreting  the  UR  wake  geometry. 

RSLTSPLT,  which  is  also  a  separate  main  program,  reads  the  con¬ 
verged  airloads,  blade-bending  moments,  bound  vorticity  strengths, 
etc.,  fix)m  the  TAPE 3  data  set  generated  toward  the  end  of  Part  2 
(see  page  48).  Using  the  plotter  routines,  RSLTSPLT  then  creates 
curves  showing  the  time  histories  as  well  as  harmonic  analyses  of 
these  converged  quantities  which  are  the  results  of  the  BLP2. 
RSLTSPLT  was  found  to  be  helpful  in  comparing  the  blade-loads 
predictions  resulting  from  a  UR  wake  geometry  with  those  from  the 
CAL  wake  geometry  and  with  experimentally  measured  blade  loads. 
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CHAPTER  4 


COMPUTED  RESULTS  AND  INTERPRETATION 


The  UR-written  wake -geometry  program  described  in  Chapter  3  was 
first  used  for  calculating  the  induced  velocity  at  a  small  number 
of  points  on  a  flat,  lifting  disc  at  an  angle  of  attack  of  2.5°. 
Figure  20  shows  the  computed  Z  -component  of  this  velocity.  The 
induced  velocity  is  upward  on  a  crescent-shape  d  area  near  the 
leading  edge  of  the  disc  and  on  an  area  immediately  behind  the 
center.  These  results  agree  well  with  the  results  of  Mangier’ s 
computations  shown  in  Figures  3  and  4. 

The  wake-geometry  program  was  then  used  to  predict  three  wake 
geometries  for  two  flight  conditions  of  the  UH-1  rotor: 

(1)  ytt  =  0.26.  (Wake  geometry  computed  from  streamlines 

emanating  from  the  flat  disc,  thus  ignoring  the  steady 
deflections  ?oc(r;>-  See  footnote,  page  61.) 

(2)  yU  =  0.26.  (Streamlines  emanating  from  the  steady 

deflected  positions  of  the  blades.) 

(3) x*  -  0.08.  (Streamlines  emanating  from  the  steady 

deflected  positions  of  the  blades.) 

The  CAL-computed  values  for  the  airloads  and  blade  flapping  moments 
are  available  (Reference  2)  for  these  flight  conditions.  (These 
were  computed  using  a  wake  geometry  input  calculated  from  a  uniform 
downwash.)  Also  available  are  experimental  measurements  for  the 
airloads  and  flapping  moments. 

The  modified  BLP2  was  then  used  to  obtain  the  blade-loads  predic¬ 
tions  resulting  from  the  UR-computed  wake  geometries.  Except  for 
the  wake-geometry  inputs,  all  other  inputs  for  these  computations 
were  identical  to  those  in  the  corresponding  CAL  runs.  Even  for 
the  wake-geometry  inputs,  the  wake  models  used  were  the  same. 


4.1  UH-1  atytt  =  0.26  (Condition  67  of  Reference  15,  Injection 

on  Flat  Disc 


Krs  =  198  feet  per  second  (forward  speed) 

OC  -  5.8  degrees  (disc  angle  of  attack) 


o<  =  2-5° 


Contours  of  zero  downwash  are  shown  by  dotted  lines.  Numbers  shown 
in  the  figure  are  X  -direction  induced  velocity  values,  positive 
downward,  normalized  on  the  foivard  velocity 

Comparison  with  Mangier' s  induced  velocity 
values  from  Figure  3 


Figure  20.  H  -Component  of  the  Induced  Velocity  on  a  Lifting 
Disc,  Using  the  UR  Program. 
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NUT  =  2  (number  of  revolutions  of  wake) 

=  0.7  (wake  advance) 

The  induced  velocity  on  the  disc,  as  computed  by  RWAKE,  showed  an 
upwash  area  near  the  forward  edge  and  just  behind  the  center  of  the 
disc.  Figure  21  shows  a  typical  instantaneous  "picture”  of  wake 
trails.  It  is  clear  that  the  trails  deviate  appreciably  from  a 
skewed  circular  helix.  Figure  22  shows  streamlines  starting  from 
90%  radius  on  the  flat  disc.  All  the  streamlines  leave  the  disc 
in  the  downward  direction,  but  those  in  the  region  /y/  >  0.8 
soon  come  under  the  influence  of  a  tip  vortex  effect,  the  ones 
farther  out  doing  so  sooner.  Streamlines  in  the  /y/<  0.8 
region  continue  downward.  The  result  is  a  separation  of  the 
streamlines  into  two  wing-tip  type  systems  emanating  from  the  tip 
areas ,  and  a  body  of  streamlines  emanating  from  the  inboard  area 
of  the  circular  wing.  The  "wing- tip”  grouping  of  streamlines 
exhibited  a  tendency  to  go  higher  (and  move  laterally)  relative 
to  the  latter  group,  which  moves  downward  comparatively  undis¬ 
turbed  by  the  presence  of  the  disc .  Similar  separation  was 
observed  in  the  water-tunnel  experiments  reported  in  Reference  16. 

The  rolling  up  of  the  tip  vortices  is  clearly  visible  in  Figures  23 
and  24,  which  show  the  progressive  distortion  of  an  initially 
circular  fluid  contour,  released  from  90%  radius  on  the  flat  disc, 
one  and  two  rotor  revolutions  after  release,  respectively. 

Figure  25  shows  a  comparison  of  the  unsteady  parts  of  the  airloads 
(time  histories),  computed  from  this  UR  wake  geometry,  with  the 
CAL-computed  airloads  and  with  experimentally  observed  values .  It 
can  be  seen  that  the  airloads  predictions  from  the  two  different 
wake  geometries  are  only  slightly  different.  The  UR  wake  airloads 
do  not  seem  to  suggest  a  decidedly  better  (or  worse)  agreement  with 
experiment. 

In  Figure  21,  a  tip- trailed  vortex  is  seen  passing  close  to  about 
the  75%  radial  station  when  the  blade  is  at  ^  =  75°.  The  induced 

velocity  at  that  station  due  to  the  nearest  tip-vortex  trail  seg¬ 
ment  was  -10.13  fps(CAL)/-3.93  fps(UR).  (The  negative  sign  indicates 
an  upward  velocity,  since  our  convention  has  positive  Z  downward. 

See  Figure  6.)  The  induced  velocity  at  this  station  due  to  the  entire 
wake  was  +9.48  f ps(CAL)/+5 . 745  fps(UR).  As  compared  to  these,  the 
normal  component  "Oig  due  to  the  control  input  was  92  0  fps  at 
that  station.  The  total  normal  component*  was  +43.74  fps(CAL) /+48 . 60 
fps(UR) . 


* ctf ,  normal  to  the  blade  chord.  See  Figure  2(a). 
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Figure  24.  Three  Views  of  a  Fluid  Contour,  Released  From  90%  Radius  on  the 
Disc,  Two  Rotor  Revolutions  After  Release  (UH-l,yU-  =  0.26). 


Conversely,  at  21.2%  radius  and  255°  azimuth,  where  significant 
differences  between  the  two  computed  airloads  are  apparent  in 
Figure  25,  the  induced  velocity  due  to  the  three  nearest  segments 
of  the  tip* trail  from  the  preceding  blade  is  found  to  be  -4.85 
fps(CAL)/-3.41  fps(UR) .  By  comparison,  the  induced  velocity  due 
to  the  entire  wake  is  -25.1  fps(CAL)/-9.918  fps(UR).  This  differ¬ 
ence  of  about  15  fps  is  probably  due  to  the  closeness  of  the 
mesh  wake,  which  curls  around  close  to  the  retreating  blade  at 
small  radii,  combined  with  the  fact  that  this  UR  wake  computation 
used  streamlines  emanating  from  the  flat  disc,  thus  placing  the 
wake  lower  than  the  CAL  computation  which  accounted  for^^  . 

The  slight  differences  between  the  airload  predictions  from  the 
two  wake  geometries  are  more  observable  at  smaller  radii.  This  is 
probably  attributable  to  the  smaller  flapping  (plunging)  velocities 
and  smaller  tangential  velocities  (and  hence  smaller  contri¬ 

bution)  at  small  radii,  which  would  make  differences  in  tne  wake- 
induced  velocity  more  noticeable  as  seen  at  the  blade  aerodynamic 

section.  Contributions  to  induced  velocity  by  specific  vortex 
filament  segments  are  given  in  Appendix  III. 

Figure  26  shows  a  comparison  of  harmonic  breakdowns  of  the  air¬ 
loads.  Figures  27  and  28  show  the  flapwise  bending  moments 
(time  histories  of  unsteady  parts  and  harmonic  components,  respec¬ 
tively)  for  the  two  computations  and  the  experimental  results. 

4,2  UH-1  at  =  0.26  (Condition  67  of  Reference  15),  Injection 
at  Steady  Deflected  Positions  of  Blades 

^5  =  188  feet  per  second  (forward  speed) 

=  5.8  degrees  (disc  angle  of  attack) 

NW  =  2  (number  of  revolutions  of  wake) 

<T  =  0.7  (wake  advance) 

The  computed  wake  geometry  for  this  case  showed  all  the  qualitative 
features  of  Figures  21  through  24.  The  tip-vortex  effect,  however, 
was  delayed  (i.e.,  became  discernible  farther  downstream)  since 
the  streamlines,  which  emanated  from  the  ’’coned”  blade  surface  this 
time,  had  to  travel  some  distance  before  they  entered  the  regime 
influenced  by  the  pressure  discontinuity  on  the  flat  disc  which  is 
the  source  of  the  tip  vortices. 

When  the  computed  geometry  was  used  as  input  to  BLP2,  the  blade 
motions  did  not  converge.  Ihe  iterative  procedure  in  Part  2 
exhibited  a  sustained-oscillations  type  of  behavior.  These  itera- 
tion-to-iteration  oscillations  in  the  generalized  blade  displacements 
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Figure?  28.  Harmonic  Analyses  of  Measured  and  Computed  Blade' 
Flapwise  Bending  Moments  (UH-1,^a  =  0.26). 


were  less  than  1.5%  for  all  structural  modes  except  the  highest 
(2nd  antisymmetric),  for  which  they  were  about  8.6%.  The  noncon- 
verged  loads,  however,  were  found  to  be  quite  close  to  the  CAL- 
computed  loads.  At  the  255°  azimuth  of  the  21.2%  radial  station, 
the  airload  prediction  agreed  better  with  the  CAL-computed  value 
than  did  the  prediction  of  Section  4.1  (see  Figure  25). 


4.3  UH-1  at ju  =  0.08  (Condition  65  of  Reference  15),  Injection 
at  Steady  Deflected  Positions  of  Blades 

^FS  =  55  .  1  feet  per  second  (forward  speed) 

ct  =  2.5  degrees  (disc  angle  of  attack) 

NW  =  6  (number  of  revolutions  of  wake) 

=  0.7  (wake  advance) 

Realizing  that  the  differences  in  the  computed  airloads  for  the 

o  .26  flight  condition  from  the  two  wake-geometry  inputs  were 
only  slight,  it  was  decided  to  carry  out  a  similar  comparison  for 
a  slower  (/*  =  0  .08)  flight  condition.  For  this  case,  the  induced 
velocity-to-forward  speed  ratio  should  be  higher  than  for  the 
faster  flight  case  (with  comparable  disc  loading) .  Airload  pre¬ 
dictions  for  this  slower  flight  condition,  therefore,  would  be 
expected  to  show  more  marked  differences  for  one  wake  geometry  as 
compared  to  the  other. 

During  the  wake-geometry  computation,  the  aownwash-to-forward 
speed  ratio  was  found  to  be  as  high  as  1.11  in  certain  regions 
toward  the  rear  of  the  disc.  This  makes  the  entire  calculation 
of  the  induced  velocity  field  (and  the  resultant  streamlines  and 
wake  geometry)  suspect,  since  the  theory  was  based  on  a  lineariza¬ 
tion  of  the  momentum  equations.  By  way  of  comparison,  the  highest 
value  of  the  downwash- to- forward  speed  ratio  was  about  0.12  for  the 
yU  =  0.26  flight  condition.  Ihus,  the  degree  of  confidence  in  the 
computed  wake  geometry  for  the^^  0.08  flight  condition  must  be 
lower  than  foTyU  =  0.26. 

As  expected,  the  upwash  area  of  the  disc  was  about  the  same  as 
before,  viz.,  a  thin  crescent  near  the  leading  edge  and  a  small 
area  behind  the  disc  center.  Many  of  the  streamlines  emanating 
from  these  regions  started  in  the  upward  direction.  (This  did  not 
happen  in  the^^  =  0.26  case  because  the  induced  upwash,  being  truly 
a  small  perturbation  to  the  fluid  velocity,  was  more  than  compen¬ 
sated  for  by  Ifcj  ’Si/iX  ,  which  is  the  normal  (downward)  component 
of  the  forward  speed.  For  the^l  =  0.26  case,  not  only  was  higher, 
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but  Oi  was  also  somewhat  larger,  5.8°,  than  the  2.5°  for  the  M  - 
0.08  flight.) 

Figures  29  and  30  show  typical  wake  trails  for  the  A  =  0.08  flight 
condition.  The  deviations  from  a  helical  wake  are  more  prominent 
in  this  case  because  of  the  higher  ratios  of  all  the  induced  veloc¬ 
ity  components  to  the  forward  speed.  Figure  31  shows  the  stream¬ 
lines  starting  from  90%  radius.  The  wing-tip  phenomenon  is  again 
evident.  (The  upward-starting  streamlines  mentioned  above  are  not 
noticeable  in  this  picture  because  the  leading  edge  crescent  with 
large  induced  upwash  was  located  outboard  of  about  92%  radius, 
and  the  upwash  region  behind  the  center  of  the  disc  extended  only 
up  to  about  50%  radius.)  This  rolling-up-of-the-tip-vortex  effect 
is  more  clearly  seen  in  Figures  32  through  35,  which  show  the 
progressive  distortion  of  an  originally  circular  fluid  contour 
released  from  90%  radius,  one  through  four  rotor  revolutions  after 
release,  respectively. 

Figures  36  and  37  show  plots  of  the  unsteady  parts  of  the  airloads 
and  bending  moments,  respectively,  for =  0.08.  It  can  be  seen 
that  the  airload  predictions  from  the  two  different  wake  geometries 
differ  significantly .  The  airload  predictions  from  the  UR  wake 
seem  to  be  in  better  agreement  with  experiment  in  the  region  around 
90°  azimuth.  While  the  rest  of  the  azimuthal  positions  show  the 
airload  histories  to  be  qualitatively  similar,  the  increase  in 
higher  harmonic  content  with  the  UR  wake  seems  clear.  This  is 
emphasized  in  the  blade  response  in  Figure  37. 

The  wake-induced  velocity  at  the  75%  radius,  when  the  blade  is  at 
if  =  75°,  is  found  to  be  14.17  fps,  and  the  total  normal  velocity 
at  that  point  is  +55.81  fps.  Thus,  for  this  slower  flight  condi¬ 
tion,  the  wake-induced  velocity  is,  as  expected,  a  more  significant 
fraction  of  the  normal  velocity  experienced  at  the  blade  section. 
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Figure  33.  Three  Views  of  a  Fluid  Contour,  Released  From  90%  Radius  on  the 
Disc,  Two  Rotor  Revolutions  After  Release  (UH-1,^  =  0.08). 
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Figure  35.  Three  Views  of  a  Fluid  Contour,  Released  From  90%  Radius  on 
Disc,  Four  Rotor  Revolutions  After  Release  (UH-lyu.  =  0.08). 


CHAPTER  5 


CONCLUSIONS  AND  POSSIBILITIES  FOR  FURTHER  WORK 


The  wake-geometry  computations  reported  in  Chapter  4  show  that  the 
flow  field  behind  a  lifting  actuator  disc  features  tip  vortices 
similar  to  those  from  a  lifting  wing. 

The  computed  airloads  for  the>#  =  0.26  flight  condition  of  the 
UH-1  teetering  rotor  seem  to  be  less  sensitive  to  wake  geometry 
than  are  those  for  the/^  =  0.08  flight  condition.  This  is  indicat¬ 
ed  by  the  slight  differences  between  the  airload  predictions  result¬ 
ing  from  two  different  wake  geometries  (CAL  rigid  helical  wake  and 
UR  distorted  wake)  for  the JLL  =0.26  condition  and  by  the  signifi¬ 
cant  differences  between  similar  predictions  for  th =  0.08 
condition.  The  accuracy  of  the  UR-computed  wake  geometry  for 
the  slower  ( yU  -  0.08)  flight  condition  is  questionable  because 
of  applying  linearized  equations  to  an  inapplicable  amplitude 
regime.  Still,  the  airload  predictions  of  the  blade-loads  program, 
even  if  this  wake  geometry  input  is  considered  "arbitrary",  are 
useful  in  providing  some  insight  as  to  the  differences  which  will 
result  from  wake-geometry  differences. 

This  stronger  sensitivity  to  wake  geometry  in  the  case  of  a 
slower  flight  condition  calls  for  some  quantitative  discussion.  As 
suggested  near  the  end  of  Chapter  4,  the  wake-induced  velocity  in 
the  case  of  the^=  0.08  condition  forms  a  more  significant  fraction 
of  the  normal  velocity  experienced  by  the  blades.  This  suggests 
an  examination  of  the  variation  of  contributions  due  to  wake  effects 
and  those  due  to  all  other  effects. 

The  faster  flight  condition  certainly  has  larger  contributions  from 
the  geometric  pitch  input  and  of  the  blade  flapping  to  the  total 
normal  velocity  at  the  blades.  Thus,  for  example,  at  the  85%  radial 
station,  the  collective  pitch  was  0.2011  rad  at  jlC  =  0.26,  as 
compared  to  0.1209  rad  at /CC  =  0.08.  The  amplitude  of  the  cyclic 
pitch  input  was  0.1195  rad  at  =  0.26,  as  compared  to  0.04  rad 
at  -  0.03.  Similarly,  the  amplitude  of  the  first  azimuthal 
harmonic  component  of  the  displacement  in  the  rigid  flapping  of  the 
teecering  blade  (normalized  on  tip  deflection)  was  0.6096  at 
y U  -  0.26,  as  compared  to  0.1830  at/M  =  0  .08.  At  the  rotor  speed 
of  33  rad/sec,  this  flapping  component  results  in  a  maximum  flapping 
velocity  of  20.1  fps  at  -  0.26,  as  compared  to  6.04  fps  at  ytt  - 
0.08.  ' 
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The  corresponding  quasi-steady  bound  vorticity  strength  compari¬ 
sons  (see  Figure  14  and  equation  (109)  )  for  the  85%  radial  station 
are  (sq  ft/Sec): 


steady  component  331.3  (yU  =  0.26)  as  compared  to  249.9  (yU  =  0.08) 
first  harmonic  158.3  (yU  =  0.26)  as  compared  to  75.34  ( jt  =  0.08) 
second  harmonic  33.2  \u  =  0.26)  as  compared  to  8.9  {yU  -  0.08) 
third  harmonic  28.91  (/U  =  0.026)  as  compared  to  13.75  (JU  =  0.08) 

These  values  show  that  in  the  case  of  the  faster  flight  condition, 
the  variations  in  the  quasi-steady  bound  vorticity  are  larger.  This 
same  point  is  illustrated  in  another  way  in  Table  IV,  which  com¬ 
pares  the  contributions  of  the  two  terms  on  the  right-hand  side  of 
equation  (109) .  The  ratio  of  the  influences  of  the  quasi-steady 
and  wake-induced  quantities  on  the  bound  vorticity  is  quite  high 
for  the  first  two  harmonic  variations  at /U  =  0.26,  but  not  so  in 
the  case  of,*  =  0.08.  For  the  3rd,  4th, and  5th  harmonics,  the 
ratio  is  smaller  at  either  value  of /t>  ,  because  there  are  no  cyclic 
variations  in  or  J>£  above  the  1st  harmonic,  and  the  blade 
flapping  displacements  are  smaller  in  harmonic  components  above  the 
2nd. 


TABLE  IV.  COMPARISON  OF  CONTRIBUTIONS  OF  THE  QUASI-STEADY 
AND  WAKE-INDUCED  QUANTITIES  TO  THE  BOUND 
VORTICITY  STRENGTH  COMPONENTS  AT  THE  85%  RADIAL 


STATION, 

UH-l,/* 

=0.26  AND/U 

=  0.08 

— 

The  Ratio 

s: 

r* 

rr— wi 

(See  Equation  (109)  ) 

Harmonic  Component 

=  0.26 

/< 

=  0.08 

Steady 

1.17 

-0.17 

1.16 

-0.16 

1st  harmonic 

1.34 

-0.34 

2.7 

-1.7 

2nd  harmonic 

1.25 

-0.25 

0.497 

0.503 

3rd  harmonic 

1.76 

-0.76 

2.94 

-1.94 

4th  harmonic 

0.684 

0.316 

0.157 

0.843 

5th  harmonic 

0.204 

0.796 

3.59 

-2.59 

As  regards  the  amount  of  wake  vorticity  remaining  close  to  the  blades 
Figures  21,  29,  and  30  show  how  much  more  rapidly  the  distance 
between  a  given  blade  and  the  wake  increases  for  the  higher  speed 
case.  The  strength  of  the  wake  vorticity,  of  course,  also  enters 
these  considerations ,  but  this  is  a  more  complex  and  less  influ¬ 
ential  variable  as  regards  the  effects  of  forward  speed. 


It  can  therefore  be  concluded  that  the  relative  insensitivity  of 
the  blade  loads  to  wake  geometry  in  the  case  of  the  faster  flight 
condition  (  -  0.26)  is  due  mainly  to  (1)  the  wake  being  blown 

away  relatively  quickly  and  (2)  a  predominance  of  the  cyclic  pitch 
input,  tangential  velocity  variation,  and  blade  flapping  displace¬ 
ments  in  the  total  time-varying  aerodynamic  environment  at  the 
blade  sections. 

It  may  be  noted  that  all  the  computations  reported  here  pertained 
to  the  UH-1  two-bladed  teetering  rotor.  It  would  probably  be 
worthwhile  to  carry  out  similar  comparative  airload  computations 
for  a  fully  articulated  rotor  such  as  the  four-bladed  Sikorsky  H-34. 
This  might  facilitate  comparisons  of  the  importance  of  blade 
displacements,  cyclic  pitch  and  tangential  velocity  variations 
and  wake-induced  effects.  Since  the  H-34  has  four  blades,  there 
will  be  more  wake  vorticity  near  the  disc,  and  hence  wake  geometry 
could  be  more  influential  than  it  seemed  in  the  jU  =0.25  case 
for  the  UH-1.  If  this  does  seem  worthwhile,  the  wake  computation 
can  be  repeated  using  the  converged  blade  loads  to  establish  the 
time-average  pressure  field  which  is  the  basic  input  to  the  UR 
wake  geometry  determination. 

If  the  UR  wake- computing  program  is  to  be  extensively  used  with 
the  CAL  BLP2 ,  considerable  saving  in  running  time  could  be 
achieved  by  writing  the  UR  wake  geometry  onto  a  temporary  direct- 
access  data  set  at  the  beginning  of  subroutine  SUB14  (see  Figure  19) 
and  reading  it  from  this  data  set  instead  of  repeatedly  reading 
it  sequentially  from  TAPE8 . 

As  alluded  to  earlier,  a  fundamental  difficulty  exists  in  tie  work 
reported  here.  Linearizing  the  equations  of  motion  on  the  basis  of 
disturbance  velocities  being  small  compared  to  forward  speed  means 
that  the  wake  computations  will  be  more  reliable,  for  a  given  tip 
speed,  at  the  higher  advance  ratios.  On  the  other  hand,  since  we 
are  using  an  assumed  steady  pressure  field  for  calculating  the  paths 
of  particles  which  are  really  unsteady,  sufficient  time  must  elapse 
to  allow  an  averaging  of  the  unsteady  effects  to  take  place.  Thus, 
the  streamline  positions--and  the  wake  geometry  obtained  therefrom-- 
will  be  most  realistic  (in  the  case  of  a  finite  number  of  blades) 
only  when  each  vorticity  element  leaving  the  blades  is  subjected  to 
the  effect  of  a  large  number  of  blade  passages.  For  a  given  number 
of  blades  and  a  given  tip  speed,  this  calls  for  a  small  advance 
ratio,*.  Put  another  way,  the  method  explored  here  will  be  most 
accurate  for  a  lightly  loaded  rotor  (i.e.,  low  thrust  per  disc 
area)  with  many  blades  at  high  forward  speed  and  high  tip  speed. 

Like  all  other  rotor  vortex  aerodynamic  theories,  the  results  will 
reflect  changes  in  vortex  strength  and  position  only  when  the  number 
of  blades  is  sufficiently  high  to  leave  vortices  close  to  succeeding 
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blades  and  yet  low  enough  to  keep  the  strength  of  the  vortices 
from  becoming  too  low.  The  scope  of  this  research  has  not 
allowed  definition  of  these  limits,  if  in  fact  they  exist.  It 
has  suggested,  however,  that  wake  details  are  important  for  tip 
speed  ratios  of  0.1  and  below,  and  are  unimportant  at  around 
=0.25  and  above. 
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APPENDIX  I 


ASSOCIATED  FUNCTIONS  OF  LEGENDRE  AND  SOME  OF  THEIR  PROPERTIES* 


1-1  Legendre  Functions  of  the  First  Kind, 


VO 

=  1 

l?co 

»  £ 

VO 

=  i-(3Z  -  O 

P3(z) 

*  *-(S£x-3') 

Z 

Particular  values: 

j 

f  Ou.  /  \  a  a** 

'!!  . . .  n  even 

( 

f 

^  •  •  •  n  odd 

s. 

*  1 

Asymptotic  values  for  |^J — *oO: 


*  Reference  24. 

**(2n-l)l[A  1«3#5*. . . .  (2n-l)  ,  with  (-!)][=  1.  (Schuster's  notation). 
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1-2  Legendre  Functions  of  the  Second  Kind,  Ogl*) 


For  Z  outside  the  segment  (-1,  -i-1) : 

r\  1  I 

«/*■>  ’T/of  jr,~ 


-l 

ae(*)‘i-(3z*-i)-ty  ~  i-* 


**(*>  $-&***);«  frr  - 


"V  ' 

For  -1  •<  Z  •«£  +1,  replaced 


foa  *  in  the  above  formulas  by  Jbg 

v*-7  / 

Particular  values  for  6+0  (6  >0)  : 

<*„<?**) 

Asymptotic  values  fov/z/+c°  '• 


aJ*) 


O! 


~ (Pn-Q  Z**7 


(Note:  from  the  above  formulas,  take  the  forms: 


Q,0?)  *  -it™7  -y 

Qi  O' 7)  =  7too~’  ’y'~1 

a2  0?)  =  -jt&v**1) ta/>  f  jr  ”  ir‘7 

ej  (‘?J  ’  -  i-(s7 "+3>  ***"■?  +T7*+  3- 


1+M 
1-Z  ’ 


and  the  values  on  the  disc  (  are: 

G,  O'o)  s  ~7 

C‘°)  =  f 

1-3  Associated  Functions  of  Legendre  of  the  First  Kind,^/»  n(z) 
For  -  7  -C  *  7  : 

P”(Z)£(-l)m(l-Z*)  m/*  (*)  (Hobson's  notation) 

dzm 

~V7^Z* 

Pi  (*)  =  -  3  *  V7^p 

pj  &•-$■  z(sz*-  o  v~n* 

Pl(*>  -  -§z(7**-3)Y7^r* 
pf(M)*3(l-z») 

p/w  -  tez(f-z'j 

For  2  outside  the  segment  (-1,  +1),  replace  (-1)*-  (1  -  E 

in  the  above  formulas  by 

Particular  values: 


_  PZ) 

P»(°) 


(-0 


n"T  (n+m-1)  // 

/iP.(2gL)' 


•  -  -(n-rr>)  ere/? 

•  .  •  (n-m)  oe/ef 
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For  e  -+0  (€>&)■ 


p:o-e)  =  c->r- 


m  /(h-sTt)/ 


PfO  +  e)  = 


{/?+m)  / 
m  ? / 


Asymptotic  values  for  /*/—  ~  : 


pf(z)  ~ 

(  n  -  /»v  / 

1-4  Associated  Functions  of  Legendre  of  the  Second  Kind,  Qn  (Z) 
For  Z  outside  the  segment  (-1,  +1): 


<27 (z)  £  (z*- 1)  m/*. 


<tmQnCz) 
dz  m 


Z-7 


i 


QjW*  Sy  JTT  -  ¥ z 


<S?i  (z)-fz  (7z*-3)  -fy  0  -  §-(S7z*-P)JiK]  -jj= 


s/(z)  -  +(*•-*)  Sy  -Jz, 

Sj(z)  =  gz(z*-,)Je*  fif  -  rgMM*  ro  *  Jf- 
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L 


For  -1  <  Z  <T  +1,  replace  the  multiplicative  factor  (Z*~ V)m^ 
in  the  above  formulas  by  (-1)®  (l-^2)*/2;  and  replace  log  -?.f  7 

by  log  JL£S  .  Z  -  1 

7-2 

Particular  values  for (6>O,*»>0): 

<(’+*)  - (-t)m 

€  * 

Asymptotic  values  fov/jr/ — 


#?(*)  ~  (-') 


rr?  (/>+/»)  /  _ ; 

C2r)  +  7)!.'  Z 


_y»  . 

Note :  From  the  above  formulas ,  Qn  take  the  forms: 


e*  (*f>  ‘ V77?1  -j£p 
*lO’r)  -Jiyi/zpto'-'t 

«’(>?)  -  /  fyvyy  - 

9+  07) '  *%<(&■?**  s)t/uJ’  <v==> 

<**0f)-3‘(t+r)ton"-L-3l7  -  ££*- 


103 


and  the  values  on  the  disc  (y  =  o)  are: 

<?j  (to)  =  -2 1 
<?;  = /  / 

<?/  -  8 


I-S  Some  Properties  of  the  Associated  Functions  of  Legendre, 

O'*)  Q"(z) 

*?(-*)  ‘(-r)m"’p?(z) 

$?(-*)  -(■Om~*’aZC*) 


Orthogonality 

1 

J  f=nm(Z)  P^(Z)dz  *0  for  nf  7,' 

o 

J  ’P"(Z)P"(Z)  yfe  -0  for  mfm' 

O  '  * 


Normalization 


/fpm(z)J  *d*  =  IrfK&fdz 

O  -r 


1  (n+m)  ? 

Pn+7  (/> -m)f 


*  d* 


(n+m)  f 
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Differentiation 


df^Cz)  = 
dz 


nz 

Z2-7 


/?+/n 

Z*-7 


PnZ(*> 


Forks'  outside  the  segment  (-1,  +1) , 


ctalcz)^ 

dz 


nZ 

Z*-7 


<*?*>  ~ 


/>+/r> 

Z*-7 


ctayov)  _ 

dy 


C/)+n>)i 
/  +  y* 
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Subroutine  VLCTS . 145 
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■jnw  etoiN  t l*  pake  r»-i e  hakm.v«ic  components  plots. 

CALL  SYMLfiL  (U»i  3  .  if  C  .2  1  ,J>Lu  TNU,  90.  ,  i.t>  I 
CALL  SYMbUL  l  ?.,  4. 8,0. ?l, HU  ««l  NUMLliC  )  ,0.  ,24) 

CALL  SYmilJL  It.-j,  9.8,  0.  14,  l  it*.  (HARMONIC  ANALYSIS)  ,0.,l9) 
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APPENDIX  III 

INDUCED  VELOCITY  AT  SELECTED  POSITIONS  ON  THE  DISC  DUE  TO 
SELECTED  SEGMENTS  OF  THE  NEAREST  TIP  TRAIL 


Table  V  shows  the  induced  velocity  due  to  nearby  filaments  of  the 
trailed  vorticity.  Five  instantaneous  positions  were  chosen  for 
points  on  the  blades  (see  Figure  38),  and  the  induced  velocities 
resulting  at  these  points  due  to  selected  short  segments  of  the 
tip  vortices  (as  positioned  by  the  two  geometries,  CAL  and  UR) 
were  computed  for  the  LL  -  0.26  flight  condition  of  UH-1. 
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Blade  at 
4i=  75° 


Blade  at 
<l>=  ZZb° 


Coordinates  (x1,  y,  z')  of  blade  stations  1,  Z,  3; 

1.  (+1.  514,  -4.434,  -0. ~44 ) 

Z.  (-3.  649,  +1-1.40,  -0.  679) 

3.  (-4. 574,  +15.86,  -0.  854) 

See  Table  V  for  coordinates  of  wake  segments  as  represented 
in  the  two  wake  geometries. 


Figure  38.  Positions  of  Selected  Blade  Stations  and  Selected 
Wake  Trail  Segments  (UH-l,yu.  =  0,26). 
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See  Figure  38. 
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